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An  airfoil  was  designed  using  the  'inverse  design*  method.  The 
airfoil  was  tailored  to  meet  the  requirements  of  high  altitude,  long 
endurance  missions.  The  cruise  Mach  number  was  0.65  so  compressibility 
effects  were  included.  In  the  procedure  used,  a  basic  airfoil  was 
generated  to  meet  the  thickness  and  moment  requirements.  It  was  then 
refined  to  ensure  desired  behaviour  of  the  boundary  layer  at  the  operat¬ 
ing  angles  of  attack.  Computer  codes  designed  by  Richard  Eppler  were 
used  for  this  study.  The  airfoil  was  anlayzed  by  using  a  viscous 
effects  analysis  program  designed  by  Shang  et  al  and  a  comparison  was 
made  with  results  obtained  through  Eppler' s  code. 


I.  INTROOXTriON 


The  United  States  Air  Force  is  interested  in  the  potential  for 
high  altitude,  long  endurance,  remotely  piloted  vehicles  to  perform 
a  variety  of  reconnaissance,  survillianoe  and  relay  missions.  The 
mission  requires  that  the  vehicle  operate  above  55,000  feet  for  at 
least  24  hours.  The  mission  is  to  be  performed  at  least  250  nautical 
miles  away  fran  hone  base.  The  aircraft  should  be  able  to  take  off 
and  land  using  conventional  runways  (Ref  11) . 

High  altitude,  long  endurance  missions  are  ideal  potential  appli 
cations  of  subcritical  airfoils.  The  wing  accounts  for  a  large  per¬ 
centage  of  total  aircraft  drag  on  such  aircraft.  This  leads  to  the 
study  of  airfoils  with  a  high  lift  to  drag  ratio.  Such  airfoils  can 
result  in  substantial  fuel  savings  when  used  on  aircraft  performing 
long  endurance  missions. 

A  substantial  amount  of  research  has  been  carried  out  in  order 
to  develop  airfoils  with  high  lift  and  low  drag  characteristics. 

But  it  has  remained  limited  to  low  speed  because  of  the  applications, 
i.e. ,  sailplanes,  gliders,  etc.  It  has  been  found  that  aircraft 
using  jet  engines  can  operate  ecncmically  at  high  subsonic  speeds. 
Thus  efforts  are  now  being  made  to  design  subcritical  airfoils  for 
high  subsonic  flight. 

Problem 

The  basic  requirement  for  this  study  was  as  follows:  To  design 
an  airfoil  with  the  following  characteristics: 

1.  A  critical  Mach  number  of  0.65. 


2.  A  nominal  thickness  ratio  of  15%. 


3.  As  high  a  design  lift  coefficient  as  possible. 

4.  A  lew  drag  coefficient  through  extensive  laminar  flow  at 
a  Reynolds  number  range  of  1.5  to  6.5  million  on  the  airfoil. 

5.  As  wide  a  drag  bucket  as  possible. 


Scope 

The  scope  of  this  study  is  well  defined  by  the  design  specifica¬ 
tions.  This  study  was  limited  to  a  theoretical  investigation  only  and 
it  did  not  include  fabrication  and  wind  tunnel  testing  of  the  airfoil . 
The  Mach  number  was  limited  to  0.65  and  the  Reynolds  number  range 
varied  frem  1.5  million  to  6.5  million.  The  requirements  of  high 
lift,  low  drag,  lav  thickness  ratio  and  a  wide  drag  bucket  indicate 
that  a  high  lift  to  drag  ratio  was  sought  in  this  design. 

General  Approach 

The  study  consisted  of  three  basic  parts  which  are  discussed 

here. 

The  first  step  was  to  use  seme  scheme  in  order  to  design  a  basic 
airfoil  shape  that  satisfied  the  design  requirements.  This  required 
the  use  of  same  oenputer  program  which  could  use  design  requirements 
as  input  to  give  the  airfoil  shape  as  output.  An  inverse  program 
designed  by  Eppler  (Ref  5)  was  used  in  this  study  to  accomplish  this 
part. 

The  second  step  involved  the  analysis  of  the  airfoil  under  viscous 
flew  conditions.  It  included  the  study  of  boundary  layer  behaviour 
in  order  to  determine  the  extent  of  laminar  and  turbulent  flow. 

Location  of  transition  point  and  study  of  flow  separation  were  also  a 


part  of  this  step.  A  viscous  flow  computer  code  designed  by  Eppler 
was  used  to  investigate  the  viscous  behaviour  of  the  airfoil. 

To  accomplish  the  two  steps  mentioned  above,  the  airfoil  was 
required  to  meet  all  the  design  requirements.  Ihis  implied  an  iterative 
solution  of  the  problem.  Various  parameters  were  changed  gradually  in 
order  to  reach  a  stage  where  all  requirements  were  fulfilled. 

The  final  step  was  to  verify  the  results  obtained  from  the  first 
tvo  steps.  The  viscous  flow  results  obtained  by  using  Eppler 's  com¬ 
puter  code  were  compared  with  another  viscous  flow  analysis  scheme.  A 
computer  program  designed  by  Shang,  Hankey  and  Dowyer  (Ref  2)  was  used 
in  this  step.  A  comparison  of  the  results  obtained  from  the  two  pro¬ 
grams  is  presented  in  the  report. 

Development 

In  this  report,  the  analysis  of  design  parameters  is  given  in 
Chapter  II.  Chapter  III  covers  the  design  procedure  and  results  fol¬ 
lowed  by  comparison  of  results  in  Chapter  IV.  Conclusions  and  recom¬ 
mendations  are  given  in  Chapter  V.  The  historical  survey  of  airfoil 
development  and  the  theory  used  in  this  study  are  included  in  the 
appendices. 


II.  ANALYSIS  OF  DESIGN  PARAMETERS 


In  the  design  process  it  is  required  first  of  all  to  decide 
upon  the  goals  that  have  to  be  achieved.  The  goals  for  this  study 
were  mentioned  in  the  introduction.  They  will  be  discussed  in  seme 

I 

detail  here. 

Mach  Number  Range.  According  to  the  specifications  for  this 
design,  the  airfoil  should  have  a  critical  Mach  nimber  of  0.65.  The 
significance  of  this  limit  can  be  understood  by  realizing  that  the 
critical  Mach  number  will  restrict  the  values  of  and  C  on  the 

Voo  p 

airfoil  surface.  Also  to  be  understood  is  the  fact  that  M  =  0.65 
falls  in  the  category  of  high  subsonic  flow  and  thus  the  compressi¬ 
bility  effects  can  not  be  neglected.  Prandtl-Glauert  relations  are 

used  to  take  care  of  these  effects  in  calculation  of  ~~  and  C  . 

Voo  P 

The  calculations  are  presented  in  Appendix  B,  where  it  can  be  seen 
that  for  a  critical  Mach  number  of  0.65  the  value  of  is  limited 

Voo 

to  1.4172  and  the  value  of  C  is  limited  to  -1.0085.  These  are  im- 

P 

portant  considerations  for  shaping  the  desired  velocity  profile.  It 
should  be  noted  that  a  lower  critical  Mach  nurber  would  result  in 
higher  values  of  and  lower  values  of  CL. 

3  Voo  p 

Reynolds  Number  Range.  The  design  specifications  lay  down  a 
Reynolds  number  range  of  1.5  to  6.5  million  for  this  airfoil.  This 
is  a  very  interesting  condition  because  the  two  extremes  of  this 
range  lie  in  two  different  design  areas.  The  1.5  million  Reynolds 
number  is  generally  considered  to  be  a  "lew"  nurber  whereas  6.5 
million  is  considered  "high" .  At  low  Reynolds  mnbers,  the  laminar 
boundary  layer  is  quite  stable.  But  an  adverse  pressure  gradient 


exists  on  the  aft  part  of  the  airfoil.  This  implies  the  existence 
of  a  turbulent  boundary  layer  in  the  aft  region.  Due  to  the  stability 
of  the  laminar  boundary  layer,  it  is  difficult  to  promote  and  con¬ 
trol  transition  at  low  Reynolds  nunbers. 

In  case  of  high  Reynolds  numbers,  the  flow  is  predominantly 
turbulent  over  the  airfoil  because  the  laminar  boundary  layer  is  very 
unstable;  thus  at  high  Reynolds  numbers  one  has  to  consider  the  ad¬ 
vantages  of  early  transition  and  existence  of  fully  developed 
turbulent  boundary  layers  while  avoiding  separation  of  flow. 

In  all  cases  the  designer  is  concerned  about  the  phenomenon 
of  boundary  layer  separation.  Separation  causes  a  great  increase 
in  the  pressure  drag  and  a  laminar  separation  also  causes  poor  initial 
conditions  for  the  turbulent  boundary  layer.  Thus  the  designer  must 
ensure  that  either  there  is  no  separation  of  flew  or  that  the  separa¬ 
tion  is  restrictred  to  the  area  close  to  the  trailing-edge  at  all  the 
flight  Reynolds  nunbers. 

Wide  Drag  Bucket.  This  design  stresses  that  the  width  of  the 
"low  drag  bucket"  should  be  as  high  as  possible  at  all  Reynolds  nunbers. 
The  physical  significance  of  the  drag  bucket  was  explained  in  Appendix 
A.  During  the  design  phase  a  wide  drag  bucket  implies  that  the  range 
of  angles  of  attack  at  which  the  flow  remains  laminar  on  the  upper 
and  the  lower  surface  be  as  large  as  possible.  This  requirement  ties 
in  with  another  phenomenon,  i.e. ,  the  thickness  to  chord  ratio  of  the 
airfoil.  It  is  explained  in  Appendix  A  that  the  thicker  airfoils  have 
a  wider  drag  bucket.  Thus  the  width  of  the  drag  bucket  is  restricted 
by  the  t/c  ratio.  The  designer  has  to  try  and  get  the  maximum  width 
consistent  with  the  t/c  ratio. 
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Thickness  to  Chord  Ratio.  A  t/c  =  15%  was  specified  for  this 
airfoil  design.  This  thickness  ratio  limits  the  width  of  the  "low 
drag  bucket"  as  explained  in  the  previous  section.  A  study  of  NACA 
laminar  flew  airfoils  indicates  that  a  drag  bucket  width  of  2°  is 
associated  with  a  t/c  of  15%.  The  designer  has  to  try  and  increase 
the  width  beyond  2°  in  order  to  shew  an  improved  design. 

High  lift.  It  is  desired  that  the  design  lift  coefficient  of 
the  airfoil  be  as  high  as  possible  consistent  with  the  other  design 
restraints.  Actually  the  application  of  this  airfoil — as  d'^scribed  in 
the  introduction— is  such  that  a  high  lift  airfoil  is  not  desired. 

The  requirement  is  for  an  airfoil  which  should  have  the  highest 
possible  V t>  ratio  at  cruise  conditions.  This  also  inplies  that 
cruise  is  the  design  condition.  Thus  the  designer  has  to  try  and  in¬ 
crease  the  ci  and  decrease  the  c<i  as  much  as  possible  in  order  to 
meet  this  requirement. 

lew  Moment.  Although  the  limits  far  the  moment  of  the  airfoil 
have  not  been  spelled  out  in  the  introduction  but  it  is  imperative  to 
have  a  lew  moment  in  order  to  avoid  the  requirements  of  a  i*rg*>  t-aii 
structure.  A  large  tail  would  result  in  unnecessary  drag  increase. 

A  lew  moment  results  from  increasing  the  laminar  length  on  the 
lower  surface  and  decreasing  the  laminar  length  on  the  upper  surface. 
But  it  has  already  been  stated  that  the  requirements  of  low  drag  imply 
a  long  laminar  length  an  the  upper  surface.  Thus  the  designer  is 
again  faced  with  conflicting  requirements  and  he  has  to  discover  the 
optimum  solution. 
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The  moment  requirement  thus  goes  into  the  decision  regarding 
the  shape  of  the  velocity  profile,  as  do  all  the  other  design  speci¬ 
fications. 

Stall  Characteristics.  It  is  important  to  design  the  airfoil  in 
such  a  manner  that  the  possibility  of  stall  within  the  operating  angle 
of  attack  range  is  eliminated.  Stall  can  occur  due  to  flow  separation 
over  the  whole  airfoil  as  a  result  of  a  leading  edge  laminar  separation 
bubble  or  due  to  the  movement  of  the  separation  point  to  the  leading 
edge.  The  latter  condition  is  avoided  through  stabilizing  the  laminar 
part  of  the  flow  and  by  controlling  the  turbulent  part  in  order  to 
restrict  separation  to  the  r<=gion  near  the  trailing-edge .  For  the 
leading  edge  laminar  bubble,  the  remedy  is  to  modify  the  velocity  pro¬ 
file  such  that  no  suction  peaks  oocur  close  to  the  leading  edge  within 
the  desired  angle  of  attack  range. 

It  can  be  seen  that  the  formulation  of  a  basic  velocity  profile 
becomes  quite  involved  due  to  all  the  design  limiatians.  The  designer 
has  to  develop  a  profile  that  will  result  in  all  the  conditions  being 
satisfied  up  to  the  desired  level. 

Basic  Velocity  Profile 

The  formulation  of  a  velocity  profile  is  the  most  important  step 
in  the  design  of  the  airfoil.  This  profile  determines  the  shape  of  the 
airfoil  and  it  also  determines  the  behaviour  of  the  boundary  layer. 

The  velocity  profile  can  be  divided  into  different  segments.  This  dis¬ 
tribution  is  explained  in  Figure  1. 

Nose  Region.  The  nose  region  includes  the  area  starting  at  the 
forward  stagnation  point  and  ending  at  about  8%  of  chord  on  the  upper 


surface.  In  this  region  the  velocity  rises  very  quickly  frcm  zero  to 
the  peak  value  at  each  angle  of  attack.  The  velocity  distribution 
here  should  be  such  that  no  leading  edge  stall  can  take  place  within 
the  operating  range  of  angles  of  attack. 

Acceleration  Region.  This  region  extends  frcm  the  end  of  the  nose 
region  up  to  the  beginning  of  the  pressure  recovery  area.  In  this 
region  it  is  desirable  to  have  favourable  or  extremely  lew  adverse 
pressure  gradients  to  maintain  laminar  flow  resulting  in  reduction 
of  skin  friction  drag. 

Pressure  Recovery  Region.  This  region  represents  the  area  where 
a  strong  adverse  pressure  gradient  exists.  It  extends  from  the  end  of 
acceleration  region  up  to  the  trailing-edge.  This  region  is  of  great 
inportance  because  the  control  of  flow  in  this  region  is  necessary  in 
order  to  have  good  lift  and  drag  characteristics.  It  is  desirable  to 
have  the  flow  transition  to  turbulent  conditions  when  entering  this 
region.  Different  people  have  came  up  with  various  schemes  to  achieve 
this.  They  include  the  introduction  of  instability  range  as  reoonroended 
by  Wortnann  (Ref  21)  and  the  turbulent  development  region  advocated 
by  Miley  (Ref  13) .  The  shape  of  the  pressure  recovery  plays  an  im¬ 
portant  part  in  determining  the  behaviour  of  the  boundary  layer.  It 
has  been  found  that  a  concave  pressure  recovery  profile  results  in 
lower  Cd  values  and  delays  separation  of  the  turbulent  boundary  layer. 

lower  Surface.  The  lewer  surface  velocity  distribution  can 


have  all  the  components  mentioned  above  or  it  can  be  only  one  com¬ 
ponent  where  the  flow  is  accelerating.  It  is  desirable  to  have  only 
accelerating  flow  on  the  lewer  surface  in  order  to  maintain  laminar  flow 


at  the  desired  angle  of  attack.  It  is  also  observed  that  a  pressure 
recovery  region  on  the  lower  surface  causes  an  undesirable  monent 
increase  cn  the  airfoil.  But  in  sane  situations  it  is  practical  to 
have  a  laminar  segment  and  then  a  pressure  recovery  region  on  the  lower 
surface.  This  is  done  to  avoid  flow  separation  on  the  lower  surface 
at  high  Reynolds  numbers. 


III.  DESIGN  PROCEDURE  AND  RESULTS 


The  design  procedure  followed  during  the  study  will  be  described 
in  this  chapter.  The  most  important  factor  in  this  procedure  is  to 
understand  the  significance  of  each  input  parameter  and  to  interpret 
correctly  the  output  of  the  program.  It  is  therefore  in  order  to 
start  by  describing  the  input  and  output  for  the  ccrnpuber  program  used 
in  the  study. 

Description  of  Input  and  Output  Parameters 

Input  Parameters.  The  input  for  the  design  part  consists  of  a 
set  of  a*  and  <j>  values  and  the  pressure  recovery  parameters  for  the 
upper  and  lower  surface.  Using  the  circle  as  the  plane  for  input, 
the  circle  is  divided  into  many  segments.  The  angle  for  each  seg¬ 
ment  is  related  to  the  chord  length  on  the  airfoil  by  the  relation 


^  (1  +  oos<j>) 

So  at  $  =  o°,  £  =  1.0  and  at  <J>  =  180,  £  =  zero.  For  each  segment  on 
the  circle  a  value  of  a*  is  specified.  This  value  signifies  the 
angle  of  attack  at  which  the  velocity  will  have  a  constant  value 
over  that  segment.  The  arc  limits — and  all  other  data  depending 
upon  arc  limits,  i.e. ,  <p.  <j>  ,  etc. — are  not  specified  in  degrees, 
rather  they  are  specified  relative  to  the  point  nutters.  The  air¬ 
foil  points  computed  are  images  of  equidistant  points  on  the  £ 
circle  (see  Figure  2) .  Any  nutter  of  points  on  the  circle  can  be 
specified  as  long  as  it  is  divisible  by  four.  So 
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Figure  2 


Segmenting  the  Interval  (0,  2ir)  for  the  Velocity  Specifi¬ 
cation  ''*(<)>)  and  the  Angle  of  Attack  Specification 
<**(<}>) 


A4>  =  27r/nc  =  360°/]^ 


where  nc  is  the  number  of  points,  =  60  was  used  in  this  study  be¬ 
cause  it  gave  an  optimum  number  of  points  on  the  airfoil  surface. 

Now  we  define 


which  is  used  to  input  arc  limits  to  the  computer.  For  example 
=  15  means  <j>^  =  90°  which  on  the  airfoil  would  indicate  the  mid 
chord  of  upper  surface.  A  set  of  v .  and  a*  values  are  presented 
in  Table  I  along  with  the  corresponding  <ju  and  x/c  values  (Ref  10) . 

Table  I 

Example  Input  Values 


4>i 

^i 

a* 

x/c 

i  =  1 

141.0° 

23.5 

8° 

0.11 

i  =  2 

165° 

27.5 

10° 

0.017 

i  =  3 

*l'*e 

0 

12° 

0.0 

i  =  4 

360° 

60 

2° 

1.0 

The  value  ^  is  specified  as  zero  but  the  exact  location  of 

the  front  stagnation  point  is  computed  by  the  program.  In  the  data 
given  above  the  circle,  and  therefore  the  airfoil,  will  have  four 
segments. 

The  second  part  of  input  consists  of  the  pressure  recovery 
specification  for  both  surfaces  of  the  airfoil.  Five  values  are 
input  for  each  surface.  The  first  number  indicates  4>g/A<{>,  i.e.. 


the  number  of  circle  segments  which  the  program  can  vary  to  return 
a  closed  airfoil.  The  second  number  specifies  the  number  of  circle 
segments  for  the  pressure  recovery,  i.e. ,  The  third  number 

specifies  the  iteration  mode  for  the  pressure  recovery.  This  nurrber 
determines  the  interpretation  of  the  fourth  and  the  fifth  number 
also.  These  iteration  modes  and  the  recovery  parameters  in  each 
mode  are  given  in  Table  II. 

Table  II 

Pressure  Recovery  Modes 

Mode  Parameter  1  Parameter  2 

0  K  y 

1  0)'  0) 

2  y  u 

In  addition  to  the  ten  numbers  specified  above,  a  design 
iteration  mode,  desired  trai ling-edge  thickness,  a  tolerance  in 
achieving  the  trail ing-edge  thickness,  and  the  integration  scheme 
have  to  be  specified.  The  detailed  description  for  these  parameters 
can  be  found  in  Ref  7. 

The  above  information  is  enough  to  solve  the  transcendental 
equation  and  find  the  airfoil  shape.  The  next  step  is  to  specify  the 
angles  of  attack  at  which  the  characteristics  of  the  airfoil  are 
desired.  The  program  can  at  this  stage  be  used  to  analyze  the 
airfoil  by  the  paneling  method.  The  input  for  boundary  layer  analysis 
consists  of  a  set  of  Reynolds  numbers  along  with  a  transition  mode. 
The  transition  modes  are  as  follows: 


1 

* 
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Mq  =  o  -  Transition  at  laminar  separation. 

Mq  =  1  or  2  -  Fixed  transition  at  given  location  on  upper  and 
lower  surface. 

Mu  =  3  -  Natural  transition  (see  equation  &-»8) 

Mu  >  3  -  Transition  with  roughness  factor  r  =  Mu~3. 

Hie  boundary  layer  calculations  are  carried  out  for  all  values 
of  angle  of  attack  specified  in  the  earlier  input. 

Program  Output.  Hie  program  output  can  be  divided  into  two 
parts,  i.e. ,  numerical  and  graphical .  The  numerical  output  consists 
of  the  airfoil  coordinates,  or  Cp  values  for  specified  angles 
of  attack,  a  detailed  printout  of  boundary  layer  development  on  the 

upper  and  lower  surfaces  for  each  set  of  Reynolds  mmber  and  angle 
of  attack,  and  a  sutnnary  of  Cp,  Qn  and  Qj  values  at  each  angle  of 
attack  and  Reynolds  number. 

The  graphical  output  consists  of  plots  of  vs  x/c,  vs 
Ci  vs  Qj,  Ci  vs  a,  Qn  vs  a,  |  turbulent  vs  x/c,  and  x/c  separation 

s» 

vs  x/c. 

A  reader  interested  in  further  details  on  program  input  and 
output  is  referred  to  Ref  7  which  contains  all  the  details  with 
examples. 

Design  Process 

The  requirements  for  this  design  imply  a  design  velocity  dis¬ 
tribution  shewn  in  Figure  3 .  It  shows  an  accelerating  flew  on 
the  whole  lower  surface .  This  is  done  to  ensure  laminar  flew  all 
along  the  lower  surface.  On  the  upper  surface  we  have  a  nose  region 
followed  by  a  long  region  of  acceleration.  This  irrplies  laminar 
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flow  over  about  80%  of  the  upper  surface  at  design  condition^.; 
Finally  on  the  upper  surface  we  have  a  pressure  recovery  region  which 
has  a  concave  shape  in  order  to  avoid  separation  of  the  turbulent 
boundary  layer.  The  pressure  recovery  parameters  are  specified 
such  that  no  laminar  separation  takes  place  prior  to  transition. 

The  nunbers  initially  chosen  for  this  type  of  profile  were: 


vx  =  8.5 

a*  =  5.0 

Recovery  area  (upper  surface) 

v2  =  27.5 

a*  =  5.0 

Acceleration  range 

V3  —  0 

a*  =  10.0 

Leading  edge  (nose  region) 

o 

VO 

II 

> 

a*  =  1 

lower  surface 

With  the  so  a  closing  contribution  X*  =  X*  =  4.0  was  used.  On  the 
upper  surface  recovery  specification  node  2  was  used  and  X  =  8.5 
was  used  with  p  =  10.0  and  ui  =  6.5.  On  the  lower  surface  no  re¬ 
covery  was  specified  so  the  X  value  was  put  equal  to  zero.  All 
other  pressure  recovery  numbers  were  realistic  on  the  lower  surface 
also.  Iteration  node  six  with  Kr  =  4  was  used  in  order  to  achieve 
a  realistic  trail ing-edge  shape. 

The  above  values  resulted  in  an  airfoil  that  was  17.38%  thick 
and  which  had  a  positive  Qrio  =  .0187.  In  order  to  bring  these  values 
to  the  specified  limits,  iterations  on  a*  values  had  to  be  carried 
cut.  The  final  values  in  the  design  process  are: 


=  8.5 

v2  =  27.5 

* 

«2 

V3  =  0 

* 

“3 

=  12.0 

V4  =  60 

aJ 

=  5.0 
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u  =  M  =  10.0 


,u)  =  a)  =  6.5 
X  =  5.0 


r# 


4 


A  =  8.5 
A*  =  X*  =  4.0 

The  result  was  an  airfoil  with  t/c  =  15.01%  ai*3  =  -.1282.  It 
can  be  noted  that  a  small  pressure  recovery  region  had  to  be 
introduced  on  the  lower  surface. 

The  shape  of  the  airfoil  and  the  velocity  distributions  for 
seme  angles  of  attack  are  shewn  in  Figure  4  .  The  airfoil  was 
given  the  number  2500. 

In  this  phase,  it  was  required  to  have  another  airfoil  with 
the  same  t/c  and  range,  which  could  be  a  competitor  when  viscous 
flow  analysis  was  carried  out.  For  the  second  design  it  was  de¬ 
cided  to  change  the  velocity  distribution  a  little  bit.  A  pressure 
recovery  region  was  introduced  on  the  lower  surface  and  it  extended 
from  mid-chord  to  the  trailing-edge.  The  upper  surface  profile 
was  to  be  the  same  as  in  airfoil  2500.  This  was  done  in  order  to 
study  the  effects  of  a  long  pressure  recovery  region  on  the  lower 
surface.  Also  it  was  intended  to  increase  the  cairber  of  the  air¬ 
foil  so  as  to  have  higher  values. 

In  this  case, again, a  lot  of  iterations  were  carried  out  to 
reach  the  final  result  and  the  final  values  are: 

Airfoil  2510  t/c  =  15.04%  =  -.1226 


vi  = 

8.5 

aj  =  5.0 

v2  = 

27.5 

a*  =  5.0 

v3  = 

0 

a*  =  10.8 

Va  ■ 

60 

a$  =  2.0 

upper  surface:  X  =  8.5,  X*  =  4.0,  y  =  10.0,  co  —  6.5  iteration  node  = 
lower  surface:  X  =  I4.5,  X*  =  4.0,  y  =  10.0,  co  =  6.5  1^=0 

Hie  airfoil  shape  with  the  velocity  profiles  is  shown  in 
Figure  5. 

Analysis  of  Airfoils.  Hie  analysis  phase  consisted  of  looing 
at  the  velocity  and  pressure  coefficient  values  on  the  airfoils  at 
different  angles  of  attack  to  determine  the  range  of  angles  of  at¬ 
tack  fit  for  cruise  at  M  =  0.65.  Next,  the  boundary  layer  behaviour 
was  investigated,  i.e.,  the  extent  of  laminar  and  turbulent  flow, 
separation  length  and  the  C1.  Ca  and  Cm  values.  All  these  values 
were  checked  for  Reynolds  number  range  of  1.5  to  6.5  million.  Hie 
detailed  results  far  each  airfoil  are  given  in  Appendix  E  and 
Appendix  F. 

Hie  section  characteristics  of  airfoil  2500  are  shewn  graph¬ 
ically  in  Figure  6  and  Figure  7.  Figure  8  and  Figure  9  contain  the 
characteristics  of  airfoil  2510.  It  can  be  seen  fran  the  figures 
mentioned  above  that  airfoil  2500  has  a  narrower  drag  bucket  width 
at  all  Reynolds  numbers  when  compared  to  airfoil  2510.  Airfoil  2500 
has  much  longer  separation  lengths  at  various  angles  of  attack  when 
ccrpared  to  airfoil  2510.  Hie  value  is  also  higher  for  airfoil 
2500. 

In  view  of  the  above  it  was  decided  to  restrict  further  work 
to  airfoil  2510  only.  On  this  airfoil  it  was  observed  that  at 
R  =  1.5  x  10®  the  airfoil  had  a  "bubble  analog"  longer  than  .03  at 
several  angles  of  attack.  This  condition  signifies  the  presence 
of  a  laminar  separation  bubble  at  the  transition  location.  In  order 


Figure  6.  SECTION  CHARACTRISTICS : AIRFOIL  2500 


SECTION  CHARACTRISTICS: AIRFOIL  2500  (cont.d) 


R=1 .5*10 


Figure  8.  SECTION  CHARACTRISTICS:AJRFOIL  2510 


Figure  9.  SECTION  CHARACTRISTICS: AIRFOIL  2510(contd) 


to  ensure  "no  laminar  separation"  the  bubble  analog  has  to  be  avoided. 

This  was  achieved  by  reducing  the  pressure  gradient  at  the  transi¬ 
tion  location.  Hie  value  of  X  was  increased  to  9.5  and  the  resulting 
profile  was  analyzed.  It  was  found  that  the  problem  of  bubble  analog 
was  completely  resolved.  It  was  also  found  that  the  section  char¬ 
acteristics  of  this  new  profile  were  satisfying  all  the  design  require¬ 
ments.  Thus  it  was  decided  to  retain  the  new  profile  as  the  final 
airfoil.  It  was  given  the  number  2700.  Hie  airfoil  shape  along  with 
velocity  profiles  at  selected  angles  of  attack  is  shown  in  Figure  10. 
Hie  section  characteristics  of  airfoil  2700  are  displayed  in  Figure 
11  and  Figure  12. 

Results.  A  summary  of  results,  obtained  by  analysis  of  airfoil 
2700,  is  presented  here.  Hie  detailed  results  are  attached  as 
Appendix  G. 

Hie  airfoil  has  a  thickness  ratio  of  15%  which  is  the  desired 
value.  It  has  a  of  -.1112  which  is  reasonably  small  and  so  it 
does  not  pose  any  problems  with  respect  to  the  tail  size. 

Hie  value  of  pressure  coefficient  remains  greater  than  -1.008  at 
angles  of  attack  ranging  from  -2.5°  to  +2.5°.  Hius  the  airfoil  can 
be  used  for  cruising  at  M  =  0.65  at  any  of  the  angles  of  attack 
within  this  range. 

Hie  drag  bucket  width  is  5°  at  a  Reynolds  number  of  1.5  million 
and  reduces  to  3°  at  a  Reynolds  number  of  6.5  million.  This  is  a 
substantial  improvement  upon  the  NACA  64  and  65  series  laminar  flow 
airfoils  that  have  a  bucket  width  of  2°  associated  with  t/x  =  15%. 

It  is  also  noted  that  the  value  of  at  the  drag  bucket  is  0.0052 
for  a  Reynolds  number  of  1.5  x  10®  and  it  reduces  to  a  value  of  .0039 


Figure fl  .  SECTION  CHAKACTHTSTICS: AIRFOIL  2700 


Figure  12.  SECTION  CHARACTRISTICS: AIRFOIL  2700 ( eon td) 


for  Rg  =  6.5  x  10®.  These  values — when  compared  to  the  NBCA  airfoils 
are  either  equal  or  lower.  Thus  the  objective  of  widening  the  drag 
bucket  has  effectively  been  met. 

The  airfoil  has  a  laminar  length  of  0.75c  on  the  upper  surface 
and  0.5c  on  the  lower  surface.  It  has  no  laminar  separation  and 
turbulent  separation  is  restricted  to  the  last  5%  of  the  chord  on 
the  upper  surface  within  the  operating  range  of  this  airfoil.  On 
the  lower  surface  no  separation  is  indicated.  These  qualities  result 
in  very  low  drag  values  and  good  lift  coefficient  values. 

Finally  one  has  to  consider  the  most  significant  factor,  i.e., 
the  I/D  ratio.  The  behaviour  of  the  1/D  value  at  different  angles  of 
attack  and  Reynolds  numbers  is  shown  in  Figures  13  -  18.  The  peak 
value  in  five  cases  is  above  100.  And  the  cruise  condition  I/D 
(a  =  2°)  is  above  or  very  close  to  100  in  all  cases.  These  are  con¬ 
sidered  good  I/D  values  for  airfoils  in  this  category. 

Figure  10  shews  the  velocity  profiles  of  airfoil  2700.  Actually 
two  sets  of  data  are  plotted  in  this  diagram,  one  from  the  design 
mode  and  the  other  from  the  analysis  mode  where  paneling  techniques 
were  used.  The  difference  between  the  two  curves  can  hardly  be 
observed.  This  implies  that  the  inverse  method  was  very  accurate 
in  converting  the  velocity  profile  into  an  airfoil  shape. 


IV.  COMPARISON  OF  RESULTS 


After  designing  and  analyzing  the  airfoil  by  using  Eppler's  com¬ 
puter  code,  it  was  decided  to  compare  results  by  running  the  same  air¬ 
foil  through  a  different  viscous  flow  analysis  scheme.  A  computer  pro¬ 
gram,  written  by  Shang,  Hankey  and  Dowyer  (Ref  2,18)  was  available.  The 
program  was  used  for  analysis  of  boundary  layers,  using  the  eddy  vis¬ 
cosity  model.  A  description  of  the  principles  used  in  this  program  is 
given  in  Appendix  H. 

The  aim  in  using  this  eddy  viscosity  model  was  to  verify  that  the 
extent  of  laminar  and  turbulent  flew  as  indicated  by  Eppler's  code 
agreed  with  results  obtained  from  this  model.  This  program  was  not 
used  for  the  calculation  of  section  characteristics  like  C^,  C^,  etc. 

The  program  input  includes  the  airfoil  coordinates,  pressure  coefficient 
at  each  airfoil  coordinate,  Reynolds  number,  altitude  and  temperature 
data,  and  some  information  about  the  desired  output.  The  output  in¬ 
cludes  the  Cf,  momentum  thickness,  displacement  thickness,  Reynolds 
numbers  based  on  momentum  and  displacement  thickness,  etc.  at  each  air¬ 
foil  station.  The  velocity  profile  of  the  boundary  layer  can  be  ob¬ 
tained  at  any  nunber  of  airfoil  stations.  Separation  of  flow  is  indi¬ 
cated  by  a  rapidly  decreasing  friction  coefficient.  But  at  separation, 
the  program  stops  further  execution,  so  the  analysis  stops  if  separa¬ 
tion  is  indicated  anywhere  on  the  airfoil.  In  this  program  the 
surface  is  approximated  by  a  flat  plate,  so  the  upper  and  the  lower 
Surfaces  have  to  be  analyzed  separately.  It  may  be  pointed  out  that 
values  of  can  be  calculated  by  vising  the  values  obtained  from 
this  program. 


Hie  analysis  was  carried  out  on  airfoil  2700.  The  test  conditions 
(input  to  the  program)  were  a  simulation  of  design  operation,  i.e. , 
altitude  of  60,000  feet  and  a  free  stream  mach  number  of  0.65.  The 
airfoil  was  numerically  tested  at  various  angle  of  attack  and  Reynolds 
nuttier  combinations. 

The  results  indicated  that  flew  separation  was  taking  place  at 
transition  on  both  surfaces  of  the  airfoil  at  all  angles  of  attack  and 
Reynolds  numbers.  On  further  analysis  it  was  found  that  adverse  pres¬ 
sure  gradient — even  of  a  lew  magnitude — caused  the  indication  of  separa¬ 
tion.  It  was,  however,  confirmed  that  there  was  no  tendency  of  leading 
edge  stall  within  the  operating  range  of  angles  of  attack.  It  was 
also  confirmed  that  the  extent  of  laminar  flow,  on  both  surfaces  of 
the  airfoil.  Hatched  very  closely  with  the  results  obtained  from 
Eppler's  program. 

It  was  concluded  that  the  eddy  viscosity  model  was  very  sensitive 
to  adverse  pressure  gradients.  The  results  of  this  analysis  reaf¬ 
firmed  the  validity  of  Eppler's  laminar  flow  analysis  method.  However, 
his  criterion  for  prediction  of  flew  separation  and  transition  to 
turbulent  flow  could  not  be  verified.  The  two  numerical  schemes  are 
based  on  different  principles.  Eppler  has  incorporated  a  lot  of 
empirical  data  in  the  solution  of  turbulent  flew  and  transition  prob¬ 
lem  (see  Appendix  D) . 

The  conflict  in  results  obtained  frem  the  two  viscious  flew  codes 
emphasizes  the  importance  of  wind  tunnel  testing.  The  actual  per¬ 
formance  of  the  airfoil  can  only  be  judged  in  the  wind  tunnel  or  by 
carrying  out  flight  tests.  These  tests  can  establish  the  validity 
of  theoretical  models. 


Table  1 


Comparison  of  Results 
ctc  =  2° 

x/c  laminar  (upper)  x/c  laminar  (lower) 


Reynolds 

Nunber 

Eppler 

Shang,  et  al 

Eppler 

Shang,  et  al 

1.5x106 

0.735 

0.731 

.505 

0.495 

2.5x10® 

0.730 

0.728 

.496 

0.490 

3.5x10® 

0.637 

0.63 

.490 

.485 

4.5x10® 

0.492 

0.485 

.486 

0.480 

5.5x10® 

0.396 

0.38 

0.482 

0.475 

6.5x10® 

0.328 

0.32 

0.474 

0.470 

V.  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

1.  This  study  was  started  with  the  intent  of  designing  an  air¬ 
foil  with  the  characteristics  specified  in  the  introduction.  This 
goal  was  achieved  and  airfoil  2700  is  the  product  of  this  study. 

2.  All  design  goals  were  met. 

3.  The  validity  of  Eppler's  design  method  was  confirmed  by  the 
use  of  paneling  techniques. 

4.  The  validity  of  Eppler's  laminar  flow  analysis  method  was 
confirmed  by  the  results  obtained  from  the  eddy  viscosity  model  de¬ 
signed  by  Shang  et  al. 

5.  The  validity  of  Eppler's  criterion  for  prediction  of  transition 
and  for  analysis  of  turbulent  boundary  layers  could  not  be  confirmed 
because  conflicting  results  were  obtained  by  using  the  computer  code 
designed  by  Shang  et  al. 

6.  The  flew  anlaysis  program  designed  by  Shang  et  al  was  found 
to  be  very  sensitive  to  adverse  pressure  gradients. 

Reoormendations 

1.  This  study  should  be  continued  in  order  to  refine  the  present 
airfoil.  This  can  be  achieved  by  using  seme  of  the  recently  developed 
theories  to  predict  the  behaviour  of  the  boundary  layer.  By  using 
these  techniques  the  airfoil  shape  may  be  modified  such  that  there  is 
no  separation  at  all  and  the  extent  of  turbulent  flew  is  optimum. 

2.  Another  interesting  topic  would  be  to  study  the  effect  of 
boundary  layer  control  on  the  characteristics  of  the  airfoil. 


The  improvement  in  C-^  and  reduction  in  Cd  may  be  substantia]  if  tech¬ 
niques  like  suction  and  blowing  are  utilized. 

3.  It  is  also  recamended  that  the  theoretical  results  be 
verified  by  building  the  airfoil  and  then  testing  it  in  a  wind 
tunnel.  That  would  be  a  test  of  the  accuracy  of  theoretical  design 
methods. 
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Appendix  A 


HISTORICAL  DEVELOPMENT 

Historical  Review 

Airfoil  design  started  when  man — in  his  desire  to  fly — copied 
the  wings  of  different  birds.  Until  early  20th  Century,  airfoil 
"design"  was  a  largely  empirical  process  drawing  inspiration  from 
birds.  In  1915  the  NACA  was  formed  and  it  started  an  experimental 
research  program.  A  huge  mass  of  data  was  accumulated  by  this  pro¬ 
gram.  The  procedure  was  to  develop  families  of  wing  sections  by  a 
systematic  experimental  approach.  During  the  experiments  a  basic 
section  was  chosen  to  start  with  and  tests  were  conducted  by  varying 
the  different  section  parameters.  Ultimately  the  desired  shape  was 
obtained  after  many  experiments.  Hie  results  of  all  experiments 
were  preserved  and  were  published  as  airfoil  section  catalogs  to  aid 
designers.  Such  work  was  also  carried  out  at  Gftttingen  in  Germany. 
During  this  period  the  theoretical  base  was  limited  to  potential 
flow  theory. 

By  the  mid  1930s  the  influence  of  viscous  scale  effects  had  been 
acknowledged,  and  boundary  layer  theory  had  been  developed.  The 
boundary  layer  concepts  were  qualitatively  incorporated  in  the  design 
of  airfoils  which  resulted  in  the  famous  NACA  6-series  "laminar 
flow"  airfoils.  These  airfoils  were  used  extensively  by  designers 
for  the  next  twenty  years.  Sane  designers  made  modifications  to 
these  designs  in  order  to  get  desired  results,  but  the  basic  material 
remained  unchanged. 

The  big  change  came  with  the  introduction  of  digital  computers 
in  the  late  1950s.  The  equations  of  potential  flew  and  viscous  flew 


theory  were  very  well  understood  at  that  time  but  they  could  not  be 
applied  fully  in  design  because  of  the  labour  involved  in  computation 
It  vas  difficult  to  solve  the  system  of  equations  resulting  from  the 
application  of  basic  equations  even  in  the  problem  of  simple  airfoils 
in  a  real  fluid.  Thus,  until  that  time,  theory  was  only  used  as  a 
guide  to  the  experimental  research  work.  In  the  wind  tunnels  it  was 
easy  to  see  "what"  happened  but  there  were  hardly  any  clues  to  "why" 
it  happened.  Thus  it  was  necessary  to  go  beyond  the  'design  by 
testing'  phase  and  to  ocme  up  with  practical,  quantitative  solutions 
to  the  equations  of  viscous  flew. 

Computers  made  it  possible  to  solve  the  oonplex  systems  of 
equations  in  reasonable  amounts  of  time.  The  comp  ter  based 
methods  ware  successfully  applied  by  Vtortmann  (Ref  2CJ  and  t’le  re¬ 
sults  were  cataloged  as  airfoils  of  the  FX-series .  Many  others  have 
since  used  the  ooirputer  not  only  to  design  airfoils  but  to  analyze 
and  test  them  under  real  flow  conditions.  Computers  led  to  the 
advent  of  what  we  call  "computational  aerodynamics "  today.  A  sur¬ 
vey  of  the  progress  of  airfoil  design  is  shown  in  Figure  19. 

(Ref  14). 

Laminar  Flow  Airfoils 

The  experimental  work  showed  that  a  laminar  boundary  layer 
could  be  maintained  on  an  airfoil  if  the  wing  surface  was  smooth 
and  a  slightly  favourable  pressure  gradient  was  maintained.  This 
was  the  theory  underlying  the  development  of  laminar  flow  airfoils. 
The  criterion  used  for  their  design  was  to  maintain  laminar  flow 
up  to  a  certain  chord  length  at  a  certain  angle  of  attack.  After 


that  point  the  pressure  was  allowed  to  increase  linearly  bade  to  its 
trailing  edge  value.  These  airfoils  attained  30  to  50  percent  smaller 
minimixn  profile  drag  than  older  airfoils  with  the  same  thickness. 

A  study  of  laminar  flow  airfoils  indicates  that  at  the  design 
angle  of  attack,  favorable  pressure  gradient  maintains  laminar  flow 

up  to  the  position  of  peak  suction  which  results  in  the  production  of 
minimum  frictional  drag.  This  is  shewn  in  the  case  of  3  NfiCA  air¬ 
foils  in  Figure  20  (Ref  15) .  The  logical  conclusion  of  this  discussion 
is  that  the  position  of  peak  suction  should  be  moved  as  far  back  as 
possible  so  as  to  have  the  minimum  possible  drag.  But  this  is  not 
true,  because  moving  the  peak  suction  back  results  in  a  steeper  pres¬ 
sure  recovery  near  the  trailing  edge  which  results  in  increased  pres¬ 
sure  drag  due  to  thicker  turbulent  boundary  layer.  It  may  even  cause 
separation.  Thus  we  conclude  that,  with  respect  to  drag,  there  is  an 
optimum  aft  limit  for  peak  suction  in  case  of  every  airfoil. 

Another  factor  to  be  considered  is  off -design  operation.  Referr¬ 
ing  to  Figure  20,  it  can  be  observed  that  at  an  angle  of  attack  of 
1.3°,  the  flow  characteristics  have  changed.  The  peak  suction  point 
rooves  forward  suddenly.  This  can  result  in  an  adverse  pressure  gradient 
causing  instability  of  the  laminar  flow  which  ultimately  leads  to 
turbulence  and  greatly  increased  frictional  drag.  If  the  behaviour 
of  laminar  flow  airfoils  is  studied  over  a  range  of  angles  of  attack, 
the  well  known  "bucket"  effect  becomes  inmediately  obvious.  There  is 
a  range  of  angles  of  attack  at  which  the  drag  has  a  minimum  value; 
above  and  below  that  range,  there  is  a  sudden  increase  in  drag.  If 
the  values  of  C,  and  C,  for  all  angles 


Figure  20. 


Velocity  Distributions  of  three  12  percent  thick 
NBCA  Airfoils  at  Two  Angles  of  Attack' (Ref.  »5) 


are  plotted,  the  curves  depict  the  "drag  bucket".  This  is  seen  in 
Figure  21  (Ref  15)  where  such  curves  (also  called  drag  polars)  for 
fair  laminar  flow  airfoils  are  drawn.  The  drag  bucket  is  the  area 
where  Qj  has  the  minimum  value.  This  phenomenon  is  absent  in  case 
of  older  airfoils  because  on  those  airfoils  the  transition  position 
changes  gradually  over  the  angle  of  attack  range. 

The  drag  bucket  concept  is  useful  to  define  a  range  of  angles 
of  attack  within  which  the  airfoil  can  operate  most  efficiently. 

But  even  within  this  range  the  point  of  maximum  efficiency  can  be 
found;  i.e. ,  the  angle  of  attack  at  which  the  C]/Ca  ratio  is  a 
maximum.  The  width  of  the  "drag  bucket"  is  closely  connected  to 
minimum  drag  and  thickness  of  the  section  in  case  of  laminar  flow 
airfoils.  It  can  be  seen  in  Figure  21  that  the  thicker  sections 
have  a  wider  drag  bucket  but  the  minimum  drag  value  is  lower  for  the 
narrow  drag  buckets.  By  going  over  such  results  for  a  large  number 
of  airfoils  the  relation  between  the  drag  bucket  width  and  the  thick¬ 
ness  becomes  very  clear.  The  width  of  the  drag  bucket  is  normally 
indicated  by  the  range  of  angle  of  attack  over  which  the  minimum 
drag  is  maintained.  For  example,  in  652~415  airfoil,  the  subscript 
'2'  indicates  that  the  difference  between  the  two  angles  of  attack — 
at  which  the  drag  starts  to  increase  frcm  its  minimum  value — is 
2  degrees.  A  study  of  such  numbers  shews  that  in  NAZA  airfoils  a 
bucket  width  of  2°  is  associated  with  a  thickness  ratio  of  15%. 

Thus  an  improvement  can  only  be  claimed  if  this  interdependence  is 
changed  such  that  a  wider  low  drag  range  is  achieved  with  the  same 
thickness  and  minimum  drag  value. 


Further  Developments 


The  study  of  laminar  flow  airfoils  was  very  useful  because  it 
led  to  the  fact  that  a  lew  drag  airfoil  is  designed  from  a  specified 
velocity  distribution.  The  velocity  distribution  is  used  to  control 
the  boundary  layer.  It  showed  the  importance  of  the  velocity  dis¬ 
tribution,  i.e.,  the  lift,  moment,  and  the  development  of  the 
boundary  layer  (drag)  are  direct  results  of  the  velocity  distribu¬ 
tion.  Also  the  velocity  distribution  defines  the  profile  of  the  air¬ 
foil. 

Analytical  techniques  have  been  developed  such  that  boundary  layer 
theory  can  be  applied  to  the  flow  about  an  airfoil  in  two  ways. 

First  the  boundary  layer  development  can  be  determined  for  a  given 
potential  flow  velocity  distribuiton.  Second,  the  potential  flow 
field,  or  at  least  some  of  its  properties,  can  be  determined  for  a 
given  boundary  layer  development.  The  second  is  known  as  the  in¬ 
verse  or  design  problem.  In  such  problems  the  potential  flow  veloc¬ 
ity  distribution  is  specified  and  the  airfoil  shape  is  computed. 

But  this  shape  has  to  be  consistent  with  the  desired  boundary  layer 
development. 

Since  the  time  of  the  development  of  NACA  laminar  flew  airfoils, 
boundary  layer  and  potential  flew  theories  have  been  steadily  im¬ 
proved.  The  flow  around  an  airfoil  can  now  be  defined  more  exactly. 
Also  the  boundary  layer  behaviour  can  be  predicted  much  more  ac¬ 
curately.  The  use  of  complex  variables  in  potential  flow  theory 
has  led  to  conformal  mapping  techniques;  also  the  use  of  digital  com¬ 
puters  has  resulted  in  fast  solutions  of  carpi  ex  boundary  layer 
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equation  systems.  Methods  are  now  available  by  which  the  boundary 
layer  transition  and  the  behaviour  of  the  turbulent  boundary  layer 
can  be  controlled. 

Control  of  Turbulent  Boundary  Layers 

It  is  understood  that  all  airfoils  have  an  adverse  pressure  on 
the  surface  close  to  the  trailing  edge.  This  adverse  pressure 
causes  transition  of  laminar  flow  to  turbulent  flow.  Thus  it  is 
imperative  to  have  a  way  of  controlling  the  development  and  behaviour 
of  turbulent  boundary  layers  so  as  to  ensure  lew  drag  penalties. 

Studies  have  shown  that  suitable  velocity  distributions  on  the  air¬ 
foil  can  result  in  turbulent  boundary  layers  which  do  not  separate 
near  the  trailing  edge  of  the  airfoil.  But  in  order  to  shape  the 
velocity  distribution  a  procedure  is  needed  through  which  separation 
of  the  boundary  layer  can  be  predicted.  Such  a  method  was  devised 
by  B.  S.  Stratford  (Ref  19)  and  it  provides  a  velocity  distribution 
which  recovers  a  given  pressure  difference  in  the  shortest  chord  length 
without  separation  of  flew.  This  velocity  distribution  is  shown  in 
Figure  22  (Ref  13)  and  it  is  compared  to  the  conventional  linear  dis¬ 
tribution  used  on  NACA  laminar  flow  airfoils.  In  Figure  22,  0  indi¬ 
cates  the  momentum  thickness  and  H  stands  for  the  boundary  layer  form 
parameter. 

Transition  Control 

On  NACA  laminar  flow  airfoils  the  velocity  distribution  is 
such  that  the  favourable  pressure  gradient  changes  suddenly  to  an 
adverse  pressure  gradient.  This  is  also  true  of  many  airfoils  with 
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Figure  22  Influence  of  the  Shape  of  the  Pressure  Recovery  Velocity 
Distribution  on  Momentum  Thickness  and  Form  Parameter 
(Ref  13) 


concave  pressure  recovery  regions.  This  abrupt  change  of  pressure 
gradient  results  in  the  creation  of  "laminar  separation  bubbles. " 

The  separation  bubble  is  a  region  where  the  laminar  flow  separates 
but  then  it  immediately  becomes  turbulent  and  reattaches  to  the  air  ¬ 
foil  surface.  The  separation  bubbles  cause  an  unnecessary  increase 
in  the  drag  by  causing  thickening  of  the  boundary  layer.  The  ef¬ 
fect  of  the  separation  bubble  can  be  seen  in  Figure  23  (Ref  15 ) . 

F.  X.  Wortman  (Ref  20)  discovered  that  it  is  possible  to  con¬ 
trol  transition  by  introducing  an  "instability  range"  between 
laminar  and  turbulent  flews.  A  slight  negative  pressure  gradient  is 
maintained  in  this  region.  This  results  in  absence  of  laminar  sep¬ 
aration  but  introduces  a  high  instability  of  the  laminar  flew.  The 
result  is  transition  frem  laminar  to  turbulent  flew  without  separa¬ 
tion. 

Wortman  has  shewn  that  the  instability  range  oembined  with  con¬ 
cave  pressure  recovery  can  enhance  the  performance  of  an  airfoil  up 
to  18  to  20  per  cent.  Comparison  of  a  NACA  laminar  flew  airfoil  and 
a  Wortman  airfoil  with  the  same  drag  bucket  width  is  shewn  in  Fig 
24  (Ref  15) .  It  can  be  seen  that  the  19.1%  thick  Wortman  airfoil 
has  18%  less  drag  than  the  18%  thick  NACA  airfoil  and  its  minimum 
lift  is  higher. 

In  the  case  of  airfoils  designed  for  very  lew  Reynolds  numbers, 

S.  J.  Miley  (Ref  13)  suggests  placing  a  turbulent  development  region 
after  the  instability  region  so  that  the  turbulent  boundary  layer 
is  fully  developed  before  facing  the  steep  adverse  pressure  gradient 
in  the  concave  recovery  region.  A  velocity  distribution  to  achieve 
this  purpose  is  shorn  in  Figure  25. 
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Figure  23  Lift-drag  Curve  of  20  percent  Thick 

Airfoil  for  two  Reynolds  Nurrbers  Illustrating 
Effect  of  Laminar  Separation  Bubble  in 
Pressure  Recovery  Region  (Ref  15) 
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Figure  25  Upper  Surface  Velocity  Distribution 
Showing  Mi  ley's  Turbulent  Development 
Region  (Ref  13) 
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ion  Procedures 


The  above  discussion  deals  with  the  properties  of  velocity  dis¬ 
tributions  over  an  airfoil.  A  very  important  objective  is  to  get 
a  suitable  airfoil  shape  from  the  given  velocity  distribution.  To  go 
slightly  further,  a  procedure  is  also  needed  to  verify  that  the  air¬ 
foil  actually  depicts  the  boundary  layer  behaviour  which  was  predicted 
by  the  velocity  distribution. 

Such  procedures  have  been  developed  and  they  apply  conformal 
mapping  techniques  for  the  design  of  an  airfoil  from  a  given  velocity 
distribution.  The  boundary  layer  analysis  is  based  on  the  current 
knowledge  available  to  the  scientist  working  on  the  problem. 

One  such  method  was  devised  by  Richard  Eppler  (Ref  5) .  The 
method  was  later  converted  into  a  ccnputer  program  which  could  be 
used  for  design  and  analysis  of  airfoils.  This  method  was  used  for 
the  present  design  project  and  it  will  be  discussed  in  detail  in  the 
next  chapter. 

Summary 

A  brief  history  of  the  development  of  airfoil  sections  has  been 
presented  followed  by  a  survey  of  the  concepts  that  led  to  present  day 
computerized  design  of  airfoils. 
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Appendix  B 


Calculation  of  Cp  and  ~r~ 

_____  .  _  V  oo 


Given:  Mcritical  =  0.65 

Unknown:  Cp  critical 
X-  critical 

Vco 

Using  the  relationship 


Cpcr  =  — —  [{ - LK^/2  ^ 

Pcr  KMcr  1  1+[(K-1)/2]M&:  “1] 

where 

Cp^j.  =  critical  pressure  coefficient 
Mcr  =  critical  Mach  nunber 
x  =  1.4  for  air 


we  get 


1.4 

Cpcx  =  (1.4) 7o.65) 2  [{i+0.2"(0?65)2}”'4  -1] 


1.2 


(Ref  12 ) 
Page  263 


Cp^  =  -1.00852 


Now  Cp  =  1  -  (^-)2 

V  2 

or  -1.00852  =  1  -  (¥-) 

*00 


so  =  1.4172 

Veo 


Thus  the  critical  Cp  and  X—  values  are: 

“  Veo 

Cp  ■  -1.00852 


,T-=  1.4172 


Appendix  C 

AIRFOIL  DESIGN  AND  ANALYSIS  MEfflCD 
It  is  understood  that  the  properties  of  an  airfoil  depend  upon 
the  potential  flow  velocity  distribution  around  it.  One  such  dis¬ 
tribution  is  shown  in  Figure  26  (Ref  13) .  This  figure  shews  the  ve¬ 
locity  distribution  on  the  upper  and  lower  surface  for  au  and  a^,  where 
a  and  a_  are  the  angles  of  attack  at  which  a  flat  plate  velocity  gradient 
is  observed  on  the  upper  and  the  lower  surface  respectively.  In  the 
insert  of  the  same  figure  it  is  shown  that  au  and  form  the  boundaries 
of  the  "drag  bucket."  To  use  the  parameters  of  this  figure  in  order 
to  design  an  airfoil,  a  suitable  inverse  solution  method  is  needed. 

That  method  should  be  able  to  ocme  up  with  an  airfoil  shape  that  will 
produce  the  specified  velocity  distributions  at  specified  angles  of  at¬ 
tack. 

Richard  Eppler  used  the  concepts  of  conformal  mapping  and  complex 
variables  to  cane  up  with  such  a  method.  He  converted  the  numerical 
scheme  into  a  computer  code.  Uiis  code  was  utilized  in  the  present 
study.  The  conputer  code  includes  a  design  program  which  yields  the  air¬ 
foil  shape  based  on  the  desired  velocity  profile.  The  second  part  of  the 
code  uses  aerodynamic  paneling  technique  to  anlayze  the  designed  airfoil. 

ihe  conformal  napping  procedure  and  Eppler' s  calculation  scheme 
will  be  discussed  in  seme  detail  in  the  following  sections.  Discussion 
of  Eppler' s  method  is  based  on  Refs  3,  5  and  13. 

Conformal  Mapping 

Plane,  irrotational,  and  incompressible  flew  can  be  represented 
by  an  analytical  function  of  the  complex  argument  Z  =  x  +  iy: 


F(Z)  =  F(x  +  iy)  =  $(x,  y)  +  iip(x,y)  (C-l) 

where  <J>  and  \p,  the  potential  and  stream  functions,  are  real  func¬ 
tions  of  x  and  y.  The  curves  <j>  =  constant  and  ip  =  constant, 
form  tvro  families  of  orthogonal  curves  in  the  xy  plane.  The  veloc¬ 
ity  components  in  the  x  and  y  direction  are  given  by 

v  _  ii  _  _ 

3x  3y  3y  3x 

The  function  F(Z)  is  called  a  ccnplex  stream  function.  From  this 
function  the  velocity  field  is  obtained  by  differentiation  in  the 
ocnplex  plane,  where 

dF _  . _ — 


=  u  -  iv  =  w  (Z) 


(C-2) 


Here  w  =  u  -  iv  is  the  conjugate  ccnplex  nurrber  to  W  =  u  +  iv,  which 
is  obtained  by  reflection  of  W  on  the  real  axis. 

For  a  circular  cylinder  of  radius  a,  approached  in  x  direction 
by  the  undisturbed  flow  velocity  the  complex  stream  function  is 


F(Z)  =  14  (Z  +  2_) 

For  a  plane  potential  vortex,  the  stream  function  is 


F  (Z)  =  In  Z 


(C-3) 


(C-4) 


where  r  is  clockwise  circulation. 

If  we  take  an  analytical  function  of  a  ccnplex  variable  and 
split  it  into  real  and  imaginary  components: 


ccrtplex  Z  plane  is  coordinated  a  point  in  the  coqplex  £  plane  that 
can  be  designated  as  the  mirror  image  of  the  point  in  Z  plane. 

When  a  point  in  the  Z  plane  moves  along  a  curve,  the  irrage  point 
moves  along  a  curve  in  the  £  plane.  The  curve  is  called  the  image 
curve  to  the  curve  in  Z  plane.  In  other  words,  the  Z  plane  is  con¬ 
formally  napped  on  to  the  £  plane  through  equation  (05) .  The 
Simplest,  and  probably  best  known,  mapping  function  is  the  Joukowsky 
mapping  function, 

,  a2 

5  =  f  (Z)  =  Z  +  ~  (06) 

It  maps  a  circle  of  radius  'a'  about  the  origin  of  the  'Z"  plane 
into  the  twice-passed  straight-line  (slit)  from  -2a  to  +2a  in  the 
£  plane. 

The  problem  of  solving  the  flow  about  a  given  body  can  be  tac¬ 
kled  by  using  conformal  napping.  Suppose  the  flew  about  a  body  of 
contour  A  in  the  Z  plane  and  its  stream  function  f  (Z)  are  known 
(usually  flow  about  a  circular  cylinder  is  assumed) .  The  flow  about 
a  body  of  contour  B  in  the  £  plane  is  to  be  found.  For  this  prob¬ 
lem  a  mapping  function 

$=f(Z)  (C-7) 

is  needed  which  can  nap  contour  A  of  Z  plane  into  contour  B  of  the  £ 
Plane.  At  the  same  time  the  known  system  of  streamlines  and  poten¬ 
tial  lines  about  the  body  A  is  being  transformed  to  the  unknown  sys¬ 
tem  of  potential  lines  and  streamlines  about  body  B.  The  velocity 
field  about  body  B  in  the  £  plane  is  found  from  the  equation. 


w«>  - 1  - 1  •  f  -  B(z>  i  - 


(C-8) 


F(Z)  and  W(Z)  are  known  fran  the  stream  function  of  the  body  A.  The 
unknown  velocity  distribution  W(£)  about  body  B  can  be  computed  from 
equation  (C-8)  if  the  mapping  function  £  =  f  (Z) ,  that  naps  body  A 
into  body  B,  is  known  (Ref  16) . 

Eppler's  Method 

The  method  applied  by  Eppler  (Ref  3)  is  based  on  conformal 
mapping.  Ihe  three  planes  involved  are  shown  in  Figure  27.  In  the 
complex  £  plane  we  start  with  infinite  parallel  flow  around  the  unit 
circle  given  by, 

£  =  £  +  in  =  r  e1*  (09) 

The  circulation  T  is  determined  such  that  the  rear  stagnation  point 
is  located  at  £  =  1.  This  flow  is  defined  by  the  complex  potential. 

F(£)  *  <j>  +  iip  =  (eiot£  +  eia  |)  -  In  £  (010) 

where  the  circulation, 

T  =  47rsin  a  (Oil) 

This  flaw  is  napped  onto  the  Z  plane  by  means  of  a  function  Z(£) , 

Z  =  x  +  iy  (012) 

If  the  inage  is  also  to  be  an  infinite  parallel  flew  with  an  angle  of 
attack  a  around  an  arbitrary  profile,  then  Z (£)  must  map  the  exterior 
of  the  £  circle  onto  the  exterior  of  the  profile  in  a  conformal 
manner.  £  =  1  must  became  the  trailing  edge  of  the  profile  and  at  the 
same  time, 

z(»)  =  co  and  (g^)*  -  real  (013) 

The  conjugate  complex  velocity  vector  W  in  the  z  plane  is 


(014) 


W  =  Ve10 


dF  _  dF/d£ 
dZ  dz7d? 


where  V  is  the  absolute  value  of  the  complex  velocity  and  9  is  the 
argument.  The  function  dF/d  is  given  by 


dF  _  -iaia-2  r  -1 
*■  2nT  C 


(015) 


Hie  problem  at  this  stage  is  to  find  the  mapping  function  Z(c) 
such  that  the  velocity  on  the  airfoil,  V,  is  a  known  function.  Hie 
function  Z (£)  is  specified  on  the  circle  £  =  e1^  and  not  on  the  air¬ 
foil,  i.e.,  not  in  terms  of  x  but  in  terms  of  V (tf> ) .  Hie  relation 
of  x  and  <p  is  defined  as 


X  =  J  (1  +  cos  ) 


(016) 


where  c  is  the  chard  length  of  the  airfoil.  This  relation  is  exact 
for  a  flat  plate  only  and  becomes  inaccurate  for  airfoils  with  thick¬ 
ness.  According  to  Eppler  (Ref  5),  if  the  function  x(<J>)  differs  too 
much  from  equation  (016)  then  an  iteration  can  be  performed  after 
specifying  a  new  function  V(<f>  (x) ) ,  but  in  practice  no  such  iteration 
has  been  found  to  be  necessary. 

If  the  velocity  V(<j>)  is  given,  then  by  taking  the  logarithm  of 
equation  (014)  we  get 

In  W  =  In  V  -  i9  =  In  In  (017) 

On  the  circle  C  =  e^.  In  ^  is  known  from  the  values  of  V($)  and 
dF  dZ 

.  Thus  the  function  which  is  the  unknown  napping  function, 
can  be  determined  using  complex  variables. 


In  considering  the  Kutta  oondiction,  and  the  singularities  in¬ 
volved  and  keeping  in  view  equation  (015) ,  dZ/dc  can  be  represented 


(018) 


in  11.(1 -i)  +  „£,  C^  +  ilWC-" 

This  form  is  introduced  because  the  power  series  does  not  diverge  on 
any  point  on  the  boundary  of  the  circle  £  =  Using  equation 
(017)  we  rewrite  the  above  equation  as 

§§  -  ln(l  -  |)]c=ei<}>  =  -  In  V(<j>)  +  ie  +  In  (^J^i* 

-  2nd  -  5i<f))  =  P(cJ>)  +  i  Q(«J>)  (019) 
dF 

Using  the  value  of  and  equation  (018)  we  get  the  equation 
for  P  (<J>),  which  is  the  real  part  of  (019), 

P<W  mJo  (am  003  +  bm  slJl  "*>  *  *  <C*20) 

The  coefficients  a^  and  bm  are  the  same  as  coefficients  of  a 
real  Fourier  series  and  they  can  be  determined  as  such.  But  in  the 
actual  computation  these  are  not  required  to  be  found.  Only  Q(<|>) 
is  required  for  the  final  solution.  New  Q(<J>)  is  equal  to 

u|0  (bm  cos  mj>  -  a^  sin  m4>)  =  ^  £2  P(t{>)  oos  d  (021) 

The  above  equation  is  used  to  determine  Q  (4>) . 

New  substituting  equation  (020)  into  equation  (019)  we  get 

ln(f)eW  -  -  in  +  *  + 

=  In  2  gPV(tf2~q~  +  iQ  +  l^(2ie1<|,/2sin  a/2)  (022) 

or 

ln(^)ei<J>  =  Pi  sin  <J>/2  |oos  (<f>/2-a)  |  ^(OV2)  (C-23) 

But  dc  =  ie2<&t>on  the  circle  so  g—  =  ie^(g|)  ei<J>  = 

■  sin  |  |oos(4>/2-a)^|yei(0f2) 

(024) 


i 
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The  above  equation  can  be  split  into  real  and  imaginary  parts  to 
yield. 


§|  -  -usin  1 1  cos  (♦/2-a)  |  cos  [|  +  Q(«J»)3  (C-25) 

^  =  -ysin  j  j  cos  (ij>/2-a)  |  sin  [j  +  Q(<J>)  ]  (C-26) 


It  can  be  seen  that  if  Q(4>)  is  known  then  the  x  and  y  coordinates 
of  the  airfoil  can  be  found. 

The  velocity  distribution  V(4>)  has  to  obey  certain  conditions. 

At  infinity  it  must  satisfy  the  condition  that  Z  (°°)  =  °°  and  (^)  =  1 . 

This  requirement  indicates  that  the  coefficient  of  c°  in  the  ex¬ 


pansion  of  equation (C-l 8) should  be  equal  to  1.  Also,  as  a  result  of 
equation (C-25) and  (C-26).  a  closed  airfoil  should  emerge.  This  con¬ 
dition  implies  that  the  coefficient  of  in  (equation  C-1R) 
should  be  zero.  Thus  from  equation  2-18  we  get 


~  =  e°+ifcto  [l+K  +  aj  +  ibi)  C1  +  .  .  .]  (C-27) 

In  order  to  satisfy  the  conditions  given  above  we  have  to  have  the 
following: 

Sq  =  0  bb  =  0 

a]_  =  1  b^  =  1 

Using  the  expansion  of  the  Fourier  series  of  equation  (C-20)  and  the 
above  values  we  get  the  following  integral  relationships: 
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(C-28) 


iraQ  =  ^  P((j>)  d<j>  =  0 


ira]_  =  £  P{<(>)  oos<f>  d$ 
2tt 

irfc>i  =  ^  P<<J>)  sin<j>  d(j> 


(C-29) 


(C~30) 


Now  it  can  be  seen  that  the  x  and  y  coordinates  of  an  airfoil  can 
be  found  using  equations (C-20) to (C-26) such  that  the  conditions  laid 
down  by  equations (C-28) to (C— 30) are  satisfied.  This  is  the  numerical 
schane  used  in  the  computer  program  where  V(<f>)  is  specified  and 
coordinates  x  and  y  are  obtained. 

Another  important  point  to  be  considered  is  that 


P<*>  '  - 


(0-31) 


in  case  we  select  a  =  a*  for  which  the  velocity  distribution  V*  (4>,a*) 
is  to  occur  on  the  airfoil.  This  function  is  independent  of  a  and 
as  a  result  we  can  use  this  function  to  obtain  the  velocity  distribu¬ 
tion  at  any  angle  of  attack,  i.e.. 


V(»,a) 


V*U,g*) 

|cos($/2-a*) 


=  f  (♦) 


(C-32) 


which  can  be  re-written  as 


V(<J>,cO  =  V*($,ct*) 


cos  (4>/2— a.) 
cos  (<J>/2-a*) 


(C-33) 


The  velocity  distribution  is  thus  specified  by  a*  and  V*(<f>,a*)  • 


As  the  right  hand  side  of  equation (C-32) is  independent  of  a  so 


it  is  not  necessary  to  select  one  a*  at  which  the  velocity  V*(4>,o*) 
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is  to  occur  on  the  airfoil.  Different  values  of  o*  can  be  selected 
for  different  segments  of  the  airfoil.  For  ary  given  segment, 
velocity  V(4>,a)  for  any  a  is  identical  to  the  specified  V*(4>,a*)  only 
if  a  =  a *(<t>) . 

It  is  necessary  to  have  a  continuous  velocity  distribution  over 
the  whole  airfoil  but  that  does  not  require  a*  (<J>)  to  be  continuous. 
Only  V*(<j>,a*)  should  compensate  for  discontinuities  of  a*U).  If 
a*(<j>)  is  given  as  piecewise  constant,  i.e.,  a*(<j>)  =  a*  for  <P 

_<  <{>i  then  at  each  discontinuity  in  a*  (<f>)  the  following  condition 
must  hold. 


which  means 


(C-36) 


In  the  computer  program  is  not  specified.  It  is  calculated 
in  the  program  such  that  equation  (C-28 ) to (C-30 ) are  satisfied  by  ensur¬ 
ing  that  condition  2-36  is  met. 


Velocity  Profile  Specification 

For  the  specification  of  velocity  distribution,  the  circle  is 
divided  into  la  arcs.  Each  arc,  i,  extends  frcm  to  Thus 
<pn  =  o  and  4>_  =  2ir.  Ihe  leading  edge  is  <J>.  „  .  For  every  arc, 

O  In  1  f  366 

V (♦,«*)  =  Vi  W($)  (0^  <  <f>  <  <^3  (C-37) 

where  W(<j>)  is  a  function  independent  of  i.  The  function  W(4>)  allows 
a  main  pressure  recovery  and  a  closure  contribution  to  be  intro¬ 
duced  on  each  surface  (Figure  ) .  It  has  the  following  form, 

1  r  2  I**1 

*<♦>  ■  t1-0-36  ]  ,g-38> 

for  the  range  0  1  <_  <p±fze  and 

a* 

«<♦>  -  [■*  {?!  I5*\r  l1-0-36!”5! W?s)  J  (C-39> 

for  the  range  4>i/£e  £  <(>  <_  2tt. 

The  property  of  the  term  in  brackets  is  such  that  if  f  (<f>)  <  0, 

{f  (<j>) }  =  o  and  if  f  (<J>)  >  o,  { f  (4>) }  =  f  (<t>) .  Thus  W(4>)  =  1  if  <f>  >■  4>s 
and  f  in  equation  (C  -38) and  if  <p  <  and  <f>  K  fs  in  equation  (C-39 ) . 

and  <p  are  always  <tt  and  and  ^  are  always  >tt. 

X  1 

Referring  to  Figure  28  it  can  be  seen  that  —  =  (l+cos<f>)  is 

e  ‘ 

used  to  give  an  impression  of  W(x)  as  given  by  equation(C-38)  and 


(C-39).  The  first  factor  in  each  equation  represents  the  main  pressure 
recovery  and  its  length  is  specified  by  ^  for  the  upper  surface  and 
(fyj  for  the  lcwer  surface.  The  total  amount  of  pressure  recovery 
and  its  shape  is  determined  by  K  and  <J>  or  K  and  ^  depending  upon 
the  upper  or  lower  surface.  The  total  amount  of  main  pressure  re¬ 
covery  on  the  upper  surface  is  given  by, 


a)  =  1  +  K  ( 


1+OOS<J>W 


ec-40) 


and  the  absolute  value  of  the  slope  of  W(x)  at  <j>  =  <j*w  is 


yK 

rd+cos^) 


(C-41) 


It  is  possible  to  specify  almost  exactly  the  shape  of  the 
pressure  recovery  by  a  careful  choice  of  K  and  y  in  equation  (C-38 ) . 
Hie  pressure  recovery  can  be  specified  in  3  different  modes: 
(0)  K  and  y  are  specified. 

(1)  id'  and  a)  are  specified. 

(2)  y  and  u  are  specified. 

These  three  modes  apply  to  the  lower  surface  also. 

The  second  factor  in  equation  2-38  and  2-39  represents  the 
closure  contribution.  It  is  normally  a  small  region  around  the 
trailing  edge.  If  %  is  >  0  then  this  term  results  in  a  convex 
velocity  decrease  as  shewn  in  Figure  28.  The  exponents  %  and 
%  are  not  specified  rather  they  are  determined  by  the  program. 


Thus  the  total  specification  of  V(<j>)  requires  only: 

(1)  The  arc  limit  <(>j_  (except 

(2)  The  values  of  at  for  every  arc. 


(3)  Pressure  recovery  specification  -<J>T.,  y,  K,  $W'  y,  K. 

w 


(4)  Hie  lengths  of  closure  contribution  as  given  by  <J>g  and  <f>g. 
This  specification  leaves  V^,  <Jk  ,  Kjj  and  free.  There 


are  I_  matching  conditions  and  equations  (C— 28)  to  (030)  must  be 

u  “ 


satisfied.  All  these  conditions  can  be  satisfied  in  closed  form 


which  leads  to  a  transcendental  equation  for  <!>•„.  The  solution  of 

1  f  X/G 


the  equation  is  obtained  by  normal  numerical  methods.  The  values  of 
K^,  and  are  easily  determined  from  the  solution  of  the  trans¬ 
cendental  equation. 

This  concludes  the  description  of  Eppler's  design  method. 


Airfoil  Analysis  Method 

The  airfoil  analysis  method  employs  aerodynamic  paneling  tech¬ 
niques.  In  paneling  the  surface  is  divided  into  a  large  nuttier  of 
very  small  panels.  Then  singularities  are  distributed  over  the 
panels.  The  flow  condition  requires  that  the  tangential  velocity 
at  the  control  point  of  each  panel  be  zero.  The  velocity  induced 
at  each  panel  is  estimated  by  use  of  available  formulae. 

The  analysis  code  uses  the  airfoil  co-ordinates  obtained  from 
the  design  method  as  the  edges  of  the  panels  (see  Figure  29) . 

The  singularities  used  are  vortices  and  they  are  distributed  para¬ 
bolical  ly  along  the  panel.  Thus,  this  program  uses  higher  order 
paneling.  The  flew  tangency  condition  is  satisfied 


Panel  between  Pn  and  Pn+1 
Shape:  Cubic  in  n(0 
Vorticity:  Parabolic 

Flow  condition:  tangential  velocity 

equals  zero  at  Pp  and  Pn+^ 


Figure  29  Panel  Method  (Ref  5) 
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at  each  airfoil  coordinate  (not  at  the  mid  point  of  the  panel  as 
is  usually  dove) .  Zero  and  90  degree  angles  of  attack  are  analyzed 
because  the  flow  for  an  arbitrary  angle  of  attack  can  be  derived 
from  these  two  solutions  by  superposition.  The  program  includes 
an  option  by  which  additional  points  can  be  splined  in  between  the 
original  airfoil  coordinates.  This  option  can  be  utilized  to  make 
the  solution  precise  in  areas  of  less  panels. 

A  brief  description  of  the  numerical  procedure  is  given  in  the 
following  pages.  This  has  been  taken  from  Ref  5 . 

The  numerical  procedure  is  based  on  finding  the  velocity  in¬ 
duced  at  a  point  (5,n)  by  the  vorticity  distribution  y (?) 
along  a  panel  which  extends  from  (  =  o  to  (  =  o  on  the  5-axis  as 
shown  in  Figure  30.  For 

YU)  -  Yj.  <1  -  £)  (C-42) 

the  velocity  components  are 
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Figure  30  Induced-velocity  Vector 
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New  for 


17 


C-45) 


(C  -46) 


(C-47) 


(C-48) 


(C-49) 


(C-50) 


Ihe  superposition  of  vorticities  in  equations (C -4 2) and (G-45) results 
in  a  linear  vorticity  distribution  with  y(?  =  o)  =  Yi  and  y(£  =  a) 

=  Yj .  The  induced  velocity  for  this  distribution  is  given  by  equa¬ 
tions  (C-43)  ,(C-44)  /C-46)  arr! (C-47)  •  In  case  of  a  parabolic  vorticity 
distribution  (eqn  C-48)  the  factor  Yp  depends  on  the  vorticities  on 
the  preceeding  and  succeeding  panels. 

If  a  panel  is  not  straight ,  it  is  divided  into  snail  subpanels, 
Hp,  in  order  to  cornpute  its  induced  velocity  at  the  point  U,n) . 

nD  is  determined  by  the  distance  —  and  — .  If  these  distances  are 

tr  a  a 
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.large a-iough  the  camber  of  the  panel  is  neglected.  If  the  distances 
are  even  greater  then  a  global  development  is  introduced  in  which 
each  panel  is  replaced  by  a  single  vortex. 

Special  higher  order  developments  are  used  for  the  velocity 
induced  at  points  on  the  panel  itself.  At  these  points,  equations 
(C-43)  XC-44),(c-46),(C-47),(C-49)and(C-50) fail.  For  the  vorticity  in 
equation  (C -4 2)  the  induced  tangential  velocity  at  the  beginning 
of  a  panel  (Figure  31)  is 

2ir  a ■  5  .1  .  TT  //".  ...  v 

—  wt  =  g  gx  +  3  g2  +  2  (C_51) 

where  g^  and  g2  are  the  slopes  shown  in  Figure  31.  For  the  vor¬ 
ticity  in  equation  (C -4  5),  the  induced  tangential  velocity  is 

Mt  =  f  91  -  |  92  (C-52) 

and  for  the  vorticity  in  equation(C-48)  it  is, 

|  Mt  -  if  91  -  |  92  <C-53> 

All  these  equations  allow  the  computation  of  the  influence  of  any 
vorticity  distribution  on  a  panel  as  specified  by  equation  (C-42 ) , 
(C-45)and(C-48)at  ary  point  P.  A  flew  condition  at  point  P  results 
in  a  linear  equation  for  the  unknown  Y.  Thus  a  system  of  rjq  linear 
equations  results,  where  nq  is  the  number  of  points  defining  the  air- 
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Appendix  D 

Laver  Analysis  Method 


The  design  and  analysis  of  airfoils  as  described  in  Appendix  C 
is  based  on  potential  flew  theory  and  it  does  not  take  into  account 
the  effects  of  viscosity.  But  in  actual  flows,  viscous  effects  are 
always  present.  Thus  the  study  of  the  behaviour  of  boundary  layers 
becomes  imperative  in  all  practical  designs. 

'Hie  boundary  layer  may  be  defined  as  a  region  of  flow,  very 
Close  to  the  surface  of  the  airfoil  (or  any  other  body  in  flow) , 
within  which  the  velocity  varies  from  zero  (in  a  body  fixed  co¬ 
ordinate  system)  to  the  full  local  stream  value  at  the  outer  edge 
of  the  boundary  layer.  The  analytical  solution  of  the  boundary 
layer  equations  will  be  discussed  in  the  following  sections  with 
particular  stress  on  the  method  used  by  Richard  Eppler  as  described 
in  Ref  4  and  5. 


Layer 


A  method  for  the  analysis  of  boundary  layers  is  described  in 
Ref  which  is  called  the  integral  method.  According  to  this  method 
if  the  tangential  velocity  component  within  the  boundary  layer  is 
U  (x,y) ,  the  length  along  the  surface  and  in  the  flow  direction  is  x, 
and  the  length  normal  to  the  surface  is  y,  then  the  potential  flew 


velocity  is 


u(x)  =  $  U(x,y> 


(D-l) 


This  displacement  thickness  is  defined  as 


m 


«l(x)  -  i  <1  -  3-> 


where  h  is  larger  than  6. 

Hie  momentum  thickness  is  defined  as 


(D-2) 


5_  (x)  =  1  (1  -  — )  —  dy 

2'  o  u  '  u  J 

e  e 


(D-3) 


And  the  energy  thickness  is  defined  as 


53<*>  -i  t1-  <S-)2) 

e  e 


(D-4) 


The  boundary  layer  shape  factors  are  defined  as 


H  -  §1 
12  62 


(D-5) 


H  =  63 
32 


(D— 6) 


Assuming  that  the  normal  velocity  at  the  surface  is  VQ  and  that 
suction  of  boundary  layer  is  negative  and  blowing  is  positive,  the 
integral  momentum  equation  can  be  written  as 


52  +  <2  +  H12>  ^  62  -  C£  +  | 


(D-7) 


Hie  integral  energy  equation  is 


6i  +  3^6.,  =  Cn  +  ^ 

3  Up  3  u  u 


(D-8) 


Hie  skin  friction  coefficient  Cf  and  the  dissipation  coefficient 
Cp  are  defined  as  follows: 


r  -  -IS 

f 


(D-9) 


C_  =  -i-  ^  t  t—  dy 
u  pu3  o  3y  J 


(D-10) 
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in  all  cases  t  is  the  shearing  stress  and  it  is  defined  as 


T  =  PV^  (D-ll) 

Nay  the  differential  equations  (D-7)  and  (D-8)  can  be  satis¬ 
fied  by  every  solution  u(x,y),  of  the  full  partial  differential  equa 
tion  for  boundary  layer  flay.  That  is  a  project  which  requires 
extensive  use  of  digital  computers  with  very  large  memories.  Ap¬ 
proximate  solutions  can,  however,  be  determined  by  allowing  velocity 
distributions  of  the  form  ~  =  f  (— ^ ,  H  (x) ) ,  where  6  is  a  thickness 
factor  and  H  is  a  shape  factor.  The  thickness  factor  pertains  to 
the  local  thickness  of  the  boundary  layer  and  the  shape  factor 
relates  to  the  local  velocity  profile.  Normally  the  shape  factor 
is  taken  into  consideration  if  it  is  required  to  find  the  velocity 
profile  within  the  boundary  layer.  So  in  normal  applications,  and 
specially  airfoil  design,  this  factor  is  neglected  and  the  distribu¬ 
tions  of  6^,  62  and  6^  are  considered. 

laminar  Boundary  Layer:  Eppler  Solution 

Using  equations  (D-10)  and  D-ll)  the  expressions  for  and 
can  be  re-written  as  follows: 


and 


Cf  = 


v  9f  62  _  e*  (H32) 
3n  p  u_  £2 

R  u_  i 


r  -  2D*  (»32) 


(D-12) 


In  this  form  it  is  clear  that  the  only  concern  is  to  solve  for 


H12^H32^'  e**H32^'  an<*  D*^H32^*  1116  Rfiynolc*s  number  is  defined  as 
R  “  arx^  the  local  Reynolds  number  based  on  momentum  thickness  is 


„  .2*2. „H_  if 

62  v  u  % 
e 


(D-14) 


for  airfoil  boundary  layers  the  reference  length  %  is  the  airfoil 
chord  length  C  and  the  reference  velocity  ue  is  the  free  stream  veloc¬ 
ity. 


Eppler  employed  solutions  of  similar  profiles  to  develop  the 
functional  relationships  for  H^2,  e*  and  D*  in  terms  of  H^-  These 
cure  presented  graphically  in  Figure  32.  The  numerical  values  are 
available  in  Ref  4  and  Ref  5.  These  relations  express  6£  and  ^3  32 
functions  of  5^  and  6^  provided  the  potential  flow  velocity  u(x) , 
suction  velocity  VQ(x) ,  and  the  Reynolds  number  are  given.  The 
shape  of  the  boundary  layer  velocity  distribution  u(y)  can  be  de¬ 
termined  from  H12  or  H32,  but  as  62  and  6^  are  involved,  is  used 
as  the  shape  factor.  The  shape  factor  for  a  Blausius  boundary  layer 
is  1.57258.  If  H32  >.  1.57258,  then  u(y)  has  no  inflection  point. 

If  H32  <_  1.57258  then  u(y)  has  an  inflection  point.  The  laminar  sep- 
aration  limit  is  set  where  ^  =  0  at  y  =  0.  At  this  condition 
H32  =  1.51509.  Also  it  can  be  verified  that  a  favorable  pressure 
gradient  and  VQ  <  o  tend  to  increase  H32  and  vice  versa. 

The  ntmerical  integration  of  equations  (D-7)  and  (D-8)  can  be 
carried  out  to  find  the  values  of  52  and  63*  Knowing  the  relationships 
of  62  and  63  with  H32,  the  value  of  H32  at  every  'x'  can  be  found. 

For  the  integration  procedure,  the  initial  values  of  S2  ar*3  63  can 


be  derived  from  closed  solutions  for  the  first  step  of  length  Ax 
(see  Ref  4,13).  Laminar  separation  is  predicted  when  =  1.51509. 
Constant  potential  flew  velocity  yields  the  constant  value  of 
H32  =  1.57258.  It  has  been  found  by  Schlichting  (Ref  16)  and  others 
that  the  approximate  solutions  obtained  from  this  method  are  quite 
close  to  the  exact  solutions. 


Turbulent 


Layer;  Eppler  Solution 


The  analysis  of  turbulent  boundary  layers  is  carried  out  in  the 
same  manner  as  laminar  boundary  layers  but  in  this  case  the  derivation 
of  expressions  for  H^*  Cf  and  is  quite  involved.  Here  selected 
velocity  distributions  are  not  used,  as  was  done  in  case  of  laminar 
flow.  Rather  some  enpirical  relations  are  employed  to  solve  this 
problem.  The  derivation  of  these  expressions  was  carried  out  by 
Wieghardt,  Ludwieg-Tillman,  and  Rotta.  The  complete  derivations  are 
available  in  Ref  4  but  only  the  final  results  are  presented  here. 


„  =  11  H32  +  15 

O  4  n 


12  48  H32  - 


c,  =  0.045716  t  (h,  0-i)  22^2]  “°*232  -1.: 


(D-15) 


(D-16) 


<T>  =  0.01  [  <H12-1) 


(D-17) 


The  above  equations  indicate  that  there  is  not  a  single  value 
of  H.j2»  for  turbulent  boundary  layers,  at  which  separation  occurs. 
According  to  this  method  it  can  only  be  said  that 

H32  >  1.58  -*■  no  separation 
H32  <  1.46  -*■  sure  separation 
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In  the  computer  program,  a  turbulent  boundary  layer  is  assured 
V7  to  have  separated  at  a  value  of  *  1.46  because  the  values 

obtained  by  the  above  procedure  are  low  under  adverse  pressure  grad¬ 
ients. 

Transition  and  Laminar  Separation  Bubbles 

It  is  important  to  predict  the  location  of  transition  from 
laminar  to  turbulent  flew  because  the  equations  are  different  for  the 
two  regimes.  Laminar  flew  is  governed  by  equations  (D-13) ,  (D—14) 
and  the  solutions  of  similar  profiles  wehreas  turbulent  flew  is 
government  by  equations  (D-15) ,  (D— 16)  and  (D-17) .  The  transition 
Criterion  used  by  Eppler  is  as  follows: 

In  R  2  >  18.4H32  -  21.74  -  0.36r  (D-18) 

where  'r'  is  the  roughness  factor.  The  idea  is  that  transition  occurs 
if  In  Rfi2  is  as  indicated  in  equation  (D-18) ,  which  was  derived 
empirically. 

The  roughness  factor  of  zero  represents  a  smooth  surface  and 
turbulence  free  initial  flow  conditions;  it  corresponds  to  natural 
transition.  Disturbances  caused  by  insects,  a  rough  surface,  or  by 
a  turbulent  free  stream  are  represented  by  r  =  4.  Roughness  factor 
of  zero  was  used  in  this  study. 

Many  times  lamj  rwr  separation  occurs  before  the  criterion  laid 
down  in  equation  (D-18)  is  satisfied.  In  such  cases  the  procedure 
is  to  switch  to  turbulent  laws  and  continue  calculation.  But,  often 
in  reed  flows,  laminar  separation  is  followed  by  turbulent 
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re-attachment  (Ref  4,5) .  The  small  area  of  laminar  separation  before 
the  turbulent  attachment  is  called  a  "laminar  separation  bubble."  Th 
adverse  effects  of  laminar  separation  bubbles  have  already  been  dis¬ 
cussed  in  Appendix  A. 

The  results  of  Eppler's  boundary  layer  method  exhibit  a  certain 
bubble  analogy.  If  the  local  Reynolds  number  is  small  enough, 
the  method  indicates  either  a  very  slowly  increasing  shape  factor 
H32  or  a  decreasing  H32  which  rapidly  reaches  the  turbulent  separa¬ 
tion  value  of  H32  =  1.46.  For  a  slowly  increasing  H32  a  bubble  ex¬ 
tending  at  least  vp  to  H32  =  1.58  should  be  expected,  ftiis  is  possi¬ 
ble  after  laminar  separation  or  in  very  steep  adverse  pressure  grad¬ 
ients  following  transitic  i.  As  the  presence  of  a  bubble  analog 
indicates  the  possibility  of  separation  in  real  flews,  so  it  is  im¬ 
portant  to  avoid  the  bubble  analog  in  oerputations. 

Section  Characteristics 

This  section  will  be  devoted  to  the  procedures  followed  by 
Eppler  (Ref  5)  in  order  to  determine  the  section  characteristics 
of  airfoils. 

The  drag  estimation  under  normal  circumstances  is  carried  out 
by  using  the  Squire-Young  formula,  i.e. , 

,ute\  5  +  H12,te 

cd  =  26  2 ,  te  (  nj  2  (I>19) 

where  'te'  refers  to  conditions  at  the  trailing  edge.  The  use  of 
this  formula  in  cases  where  the  f lw  at  the  tr ailing-edge  is  near 
separation,  results  in  very  high  c^  values.  Actually  (D-19)  was  not 


derived  to  cater  for  such  conditions.  Empirical  modifications  were 
carried  out  to  cover  such  situations  and  the  result  is: 


where 


c„  =  2S2,  ts  (HtS)  5  +  H*lS.te 
00 


(D-20) 


H*12,te  “  H12,te  for  H12,te  -  2,5 

=  2.5  for  H. -  .  >  2.5 

12, te  — 


(D-21) 


Hie  above  relation  is  applied  to  the  trailing-edge  for  upper  and 
lower  surfaces  separately. 

In  case  of  separation,  the  momentum  thickness  62  is  corrected 


as  follows: 


52,te  *  «2,sep 


(D-22) 


where  u  ___  is  the  local  velocity  at  which  separation  occurs.  As 


H,~  =  1.46  for  turbulent  separation  so  H..  =  2.803  and  as  such 

iz,sep 

cd  frcm  equation  (D-19)  will  be 


c ,  =  c,  (H_SSE) 
d  d,sep  u  te 


Hie  above  expression  shews  that  a  penalty  in  the  form  of  drag  in¬ 
crease  results  from  separation. 

Lift  and  pitching  moment  coefficients  are  basically  determined 
by  potential  flew  methods  but  seme  viscous  effect  corrections  are 
applied  to  them. 

For  no  separation  the  lift  curve  slope  is  kept  at  2n.  Hiis  means 
that  the  potential  flow  thickness  effects  are  offset  by  the  boundary 
layer  displacement  effects.  So  if  aQ  is  the  zero  lift  angle  as  shown 
in  Figure  33,  then  the  lift  coefficient  with  no  separation  will  be 
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(D-23) 


CJlns 


=  2n(ac  +  aQ) 


where  a  is  the  angle  of  attack  with  respect  to  the  chord  line. 

In  case  a  separation  is  predicted  on  the  upper  surface  and  its 

length  is  S  then  an  angle  of  attack  correction  is  applied, 
sep 

Aa  -  -  1/2Ssse  (5us  +  ac)  (D-24) 

where  6^  is  the  slope  of  the  upper  surface  near  the  tr ailing-edge 
(see  Figure  33) . 

The  lift  coefficient  correction  will  be 


Ac^  =  2rrAa  =  -tt  (6^  +  «c)  (D-25) 

Similar  corrections  are  applied  if  separation  is  indicated  on 
the  lower  surface.  A  positive  Ac^  on  the  upper  surface  or  a  negative 
Ac^  for  the  lower  surface  are  not  allowed,  because  they  are  un¬ 
realistic. 

Summary 

The  basic  boundary  layer  concepts  were  briefly  covered  followed 
by  derivation  of  boundary  layer  solution  equations.  In  the  solu¬ 
tions,  special  emphasis  was  laid  on  Eppler 's  solution  methods  as 
they  are  of  primary  interest  in  this  project.  The  chapter  was  con¬ 
cluded  with  the  estimation  procedures  used  to  determine  section 
characteristics  like  c„,  c,  and  c  . 


Appendix  E 

Detailed  Results:  Airfoil  2500 

The  detailed  results  for  airfoil  2500  are  presented  in  this  ap¬ 
pendix.  On  the  first  two  pages,  the  x  and  y  coordiantes  of  the  air¬ 
foil  are  given  under  aolumns  marked  x  and  y.  The  Cp  values  at  each 
coordinate  are  listed  in  the  subsequent  columns.  Every  column  indi¬ 
cates  Cp  values  for  a  different  angle  of  attack.  All  angles  of 
attack  are  with  respect  to  the  chord  line.  The  last  line  on  page  9<f 

indicates  the  value  of  zero  lift  angle  of  attack  and  C  . 

mQ 

The  last  four  pages  of  the  appendix  contain  a  summary  of  the 
results  obtained  by  viscious  flew  analysis.  For  every  Reynolds 
number  (1.5  to  6.5  million)  and  angle  of  attack  (-5°  to  +11°), 
the  amount  of  turbulent  flow,  amount  of  separation,  and  the  value  of 
Cd  is  given  for  the  upper  and  lower  surfaces,  separately.  The  values 
of  turbulent  flow  and  separation  are  given  as  fractions  of  the  chord 
length,  i.e.  s.sep  =  .03  implies  that  there  is  separated  flow  on 
3%  of  chord  length. 

The  values  of  C^,  and  for  each  set  of  Reynolds  number  and 
angle  of  attack  are  also  listed  in  these  data. 
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Appendix  F 

Detailed  Results:  Airfoil  2510 

The  detailed  results  for  airfoil  2510  are  presented  in  this  ap¬ 
pendix.  On  the  first  two  pages,  the  x  and  y  coordinates  of  the  airfoil 
are  given  under  oolirms  narked  x  and  y.  The  Cp  value  at  each  co¬ 
ordinate  are  listed  in  the  subsequent  columns.  Every  column  indi¬ 
cates  Cp  values  for  a  different  angle  of  attack.  All  angles  of 
attack  are  with  respect  to  the  chord  line.  The  last  line  on  page  10} 
indicates  the  value  of  zero  lift  angle  of  attack  and  Cfc^. 

The  last  four  pages  of  the  appendix  contain  a  surrmary  of  the 
results  obtained  by  viscous  flow  analysis .  For  every  Reynolds 
nunber  (1.5  to  6.5  million)  and  angle  of  attack  (-5°  to  _11°),  the 
amount  of  turbulent  flow,  amount  of  separation,  and  the  value  of 
is  given  for  the  upper  and  lower  surfaces,  separately.  The  values 
of  turbulent  flow  and  separation  are  given  as  fractions  of  the  chord 
length,  i.e.  s.sep  -  .03  implies  that  there  is  separated  flew  on  3% 
of  chord  length. 

The  values  of  C^,  Cj  and  for  each  set  of  Reynolds  nunber  and 
angle  of  attack  are  also  listed  in  these  data. 
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Appendix  G 

Detailed  Results:  Airfoil  2700 

Hie  detailed  results  for  airfoil  2700  are  presented  in  this  ap¬ 
pendix.  On  the  first  two  pages,  the  x  and  y  coordinates  of  the  airfoil 
are  given  under  columns  marked  x  and  y.  Hie  values  at  each  co¬ 
ordinate  are  listed  in  the  subsequent  columns.  Every  column  indi¬ 
cates  Cp  values  for  a  different  angle  of  attack.  All  angles  of  at¬ 
tack  are  with  respect  to  chord  line.  Hie  last  line  on  page  100 
indicates  the  value  of  zero  lift  angle  of  attack  and  Qi^. 

Hie  last  four  pages  of  the  appendix  contain  a  sumary  of  the 
results  obtained  by  viscous  flow  analysis.  For  every  Reynolds 
number  (1.5  to  6.5  million)  and  angle  of  attack  (-5°  to  +11°),  the 
amount  of  turbulent  flow,  amount  of  separation,  and  the  value  of 
is  given  for  the  upper  and  lower  surfaces,  separately.  Hie  values  of 
turbulent  flow  and  separation  are  given  as  fractions  of  the  chord 
length,  i.e.  s.sep  «  .03  implies  that  there  is  separated  flow  over 
3%  of  chord  length. 

Hie  values  of  C^,  and  Cn^  for  each  set  of  Reynolds  number  and 
angle  of  attack  are  also  listed  in  these  data. 
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Appendix  H 


Brief  Description  of  Boundary  Layer  Analysis 
Method  Designed  by  Sharxj  et  al 

A  ccnputer  code,  using  the  eddy  viscosity  model  for  analysis 
of  boundary  layers,  was  employed  for  verification  of  results 
achieved  by  using  Eppler’s  boundary  layer  analysis  method.  A  brief 
description  of  the  method  used  in  solving  the  boundary  problem  is 
given  here. 

The  method  is  based  on  the  solution  of  three  basic  equations 
i.e.  continuity,  momentum,  and  energy  equations.  They  are: 

Continuity: 


(ri  pu)  +  §-  [ri  p(V 


P 


=  0 


(H-l) 


Momentum: 


r  9u  ...  o  'V'  *  3Ui 

p[u5J+  "  + V-’  5y  = 


dP  .  1  3  r_i  ,  3u 


®  +  57Ir:>(wi7+pu,v,,)  (H_2) 


Energy: 


p[u  ^  (CJT)  +  (V  +  ^-)  4  (C„T)]  =  u&  +  ^  £[ri  ^(CJT)] 


3u 

57  VV 


sy  p 


^(57)2  +  4I?[i:3  <-cPpv'T'»  -pu’ 


dx  rj  9y  Cp  3y  p 
(H-3) 


y. 

3y 


These  equations  were  modified  by  using  the  following  factors 

V  =  y  +  P—V— 

P 

c  _  .  u'V' 

e  ‘  p  5u75y 

V'T' 

"  “S P  5r75y 


(H-4) 
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Also  included  were  the  following  two  laws: 
Perfect  gas  law: 

p  =  Cp  (£^)p  T 
Sutherland  viscosity  law; 


(H-5) 


y  ,T  ,1.5  Te+S 

ye  V 


ff|  (air  only) 


(H-6) 


The  system  thus  has  two  algebraic  expressions  and  three  non¬ 
linear  partial  differential  equations.  But  it  still  had  a  singular¬ 
ity  at  x  =  0  (leading  edge) .  This  was  avoided  by  using  a  transforma¬ 
tion  of  variables.  Ihe  transformation  of  Prcbstein-Elliott  and 

Levy-Lees  is  used.  It  can  be  written  as: 

x 


1  - 

5  (x)  -  I  pe  Ue/ye  rp  dx 

(H-7) 

n  (X,y>  -  EsSsS^  hi  P  ^ 

2?  o  Pe 

(H-8) 

■ 

(sl’y'PeSPer^l^),,*  (fjj)  ($£ 

(H-9) 

,9)  .  Peuerljtj  .  p. .  3. 

Vx  H  W  Vs 

(H-10) 

■4 

- 

f 

3i# 

n 

(H-ll) 

i-fc 

(H-12) 

1 

(H-13) 

l  _  py 

PeVe 

(H-14) 

■ 

4 

ue2 
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CpTe 

(H-15) 
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g  _  ii 

ue  d? 


(H-16) 


Continuity; 

|£+85  H+F=0  (H-17) 

Momentum: 

8£  F  H +  v  If  -  sf  (t2j  1 1  lr,+  6  (f2  -  J  - 0  «-«> 

Energy; 

8£F  If  ■ ♦  < -  £  <‘23  k  « Is*  -  tt2j*<!?>2  * 0  «-19» 

where  e  =  1  +  -r  and  e  =  1  +  ^-  T 
V  U  Prt 

These  equations  are  cast  into  a  finite  difference  scheme  and  are 
used  to  study  the  behaviour  of  the  boundary  layer.  The  transformation 
frcm  the  real  (x,y)  plane  to  the  (5,n)  plane  cast  the  boundary  layer 
into  a  rectangular  grid  of  nodes.  The  value  of  F,  £  and  V  was  com¬ 
puted  at  each  node  in  the  grid.  This  procedure  was  carried  out 
along  the  surface  of  the  test  object  (e.g.  airfoil) .  Based  on  the 

N 

results  obtained  at  each  node  the  following  quantities  were  confuted: 
Boundary  layer  mcmentun  thickness 
Boundary  layer  displacement  thickness 
Local  skin  friction  co-efficient 
Local  Mach  nunber 
Local  Reynolds  number 

Reynolds  nunber  based  on  momentum  thickness 
Reynolds  nunber  based  on  displacement  thickness 
Ratio  of  local  and  free  stream  density  of  air 
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Ratio  of  local  and  free  stream  temperature 

Stanton  number  (for  heat  transfer) 

These  parameters  are  enough  to  enable  a  thorough  analysis  of 
the  boundary  layer  at  any  point  on  the  surface.  The  computer  code 
also  has  the  capability  of  providing  the  velocity  and  temperature 
profile  within  the  boundary  layer  at  any  number  of  points. 

The  input  to  the  program  includes  the  specification  of  free 
stream  Mach  number,  Reynolds  number  and  temperature  in  addition  to 
the  ratio  of  specific  heats,  the  laminar  and  turbulent  prandtl 
number,  the  number  and  location  of  points  (on  the  surface)  on 
which  data  is  needed,  and  the  specification  the  Cp  value  at  each 
point  on  the  surface. 

Readers  interested  in  further  details  of  this  program  are 
referred  to  Ref  18  and  Ref  2. 
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